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OBJECTIVES 


This  project  involved  theoretical  calculations  of  microscopic  atomic  and  molecular 
processes  important  in  the  generation  of  H’  in  a  hydrogen  plasma.  Efficient  production  of 
H"  is  essential  for  many  practical  applications  such  as  particle  beams  and  heating  fusion 
plasmas.  The  calculations  performed  provided  detailed  information  about  the  effect  of  the 
internal  energy  of  the  reactants.  Such  data  are  necessary  to  determine  the  optimal 
conditions  for  operating  the  plasma.  In  more  general  terms,  the  theoretical  analysis  and 
close  interaction  with  experiments  led  to  a  better  understanding  of  the  basic  physics  of  the 
plasma. 

The  specific  objectives  of  the  project  were  (l''  to  calculate  cross  sections  for 
dissociative  recombination  of  electrons  with  H2  for  various  values  of  the  vibrational 
quantum  number  v,  (2)  to  provide  theoretical  analysis  for  ongoing  experiments  using 
multiphoton  ionization  to  investigate  the  spectroscopy  of  excited  states  of  molecular 
hydrogen,  and  (3)  to  calculate  cross  sections  for  dissociative  attachment  of  electrons  with 
vibrationally  excited  H2  molecules.  The  following  sections  report  the  accomplishments  of 
all  these  investigations. 
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RESEARCH  ACCOMPLISHMENTS 


ELECTRON  RECOMBINATION  WITH  Hj 

Dissociative  recombination  involves  the  capture  of  an  electron  by  a  molecular  ion, 
formation  of  a  compound  state,  and  subsequent  dissociation  into  fragments.  The  process 
may  be  written  schematically  as 

e  +  H2  ->  H2*  ->  H*  +  H 

Initial  calculations  performed  for  e  +  H2  illuminated  the  mechanisms  that  come  into  play  as 
the  process  unfolds.  Our  calculations  treated  in  detail  the  so-called  indirect  mechanism. 

The  calculations1  were  published  in  1988,  and  a  copy  of  the  work  is  included  as 
Appendix  A. 

We  have  an  ongoing  collaboration  with  Professor  Brian  Mitchell  at  the  University 
of  Western  Ontario,  and  in  1990  we  compared  our  calculations  with  the  most  recent 
experiments  of  his  group.  The  new  data  show  dramatic  oscillations  in  the  cross  sections  as 
a  function  of  the  electron  energy.  These  results  are  of  sufficient  accuracy  to  associate 
particular  features  of  the  cross  sections  with  specific  compound  states  of  the  neutral, 
excited  hydrogen  molecule.  Additional  theoretical  analysis  was  performed  to  allow 
comparison  of  our  calculations  with  experiments  at  this  level  of  accuracy. 

Our  new  work  clearly  shows  that  both  so  and  da  Rydberg  series  of  H  2  contribute 
to  the  oscillations  in  the  cross  section,  and  we  were  able  to  make  tentative  assignments  of 
the  quantum  numbers  of  these  states  that  contribute  to  the  observed  resonances.  A  joint 
manuscript  with  Professor  Mitchell's  group  has  been  submitted  to  Physical  Review  Letters, 
and  a  copy  is  included  as  Appendix  B.  The  excellent  agreement  achieved  between  theory 
and  experiment  is  shown  in  Figure  1 . 

MULTIPHOTON  IONIZATION  OF  H2 

Two  experiments  performed  in  the  Molecular  Physics  Laboratory  at  SRI  are  related 
to  the  detection  of  vibrationally  excited  hydrogen  by  multiphoton  ionization.  This  detection 
is  an  important  diagnostic  for  hydrogen  plasma  sources  of  H\  Theoretical  calculations 
have  made  key  contributions  to  the  analysis  of  these  experiments. 
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RA-5549-26 

Comparison  of  theory  and  experiment  for  dissociative  recombination  of 
electrons  with  Hi. 


The  first  set  of  experiments  was  performed  in  Dr.  William  K.  Bischel's  research 
group.  Absolute  cross  sections  were  measured  for  two-photon  excitation  of  molecular 
hydrogen  to  the  EF  state  and  for  subsequent  ionization  of  this  state.  Theoretical 
calculations  were  performed  to  clarify  the  mechanism  of  the  ionization  of  the  EF  state. 
These  calculations  were  recently  published,2  and  a  copy  is  attached  as  Appendix  C.  The 
main  conclusion  was  that  an  autoionizing  state  provides  an  important  indirect  pathway  for 
photoionization.  The  additional  photoionization  signal  caused  by  the  indirect  ionization  is 
frequency-dependent.  Understanding  these  effects  is  necessary  for  plasma  diagnostics, 
because  the  strength  of  the  photoionization  signal  as  a  function  of  frequency  must  be 
correlated  with  the  population  of  various  vibrationally  excited  states  of  H2. 

The  second  set  of  experiments  was  performed  by  Dr.  Ravinder  Kachru,  Dr. 
Hanspeter  Helm,  and  their  associates.  The  same  multiphoton  ionization  process  of 
molecular  hydrogen  was  studied.  The  energy  distribution  of  the  photoelectrons  was 
measured,  and  calculations  of  this  energy  distribution  were  compared  with  the  experimental 
data.  Good  agreement  was  achieved,  leading  to  confidence  in  our  understanding  of  the 
physics  of  the  process.  A  joint  experimental  and  theoretical  paper  was  published3 
describing  the  work,  and  a  copy  of  this  paper  is  attached  as  Appendix  D. 

ELECTRON  ATTACHMENT  TO  H2 

The  process  of  electron  attachment  to  h2  may  be  written 

e  +  H2  ->  H'  +  H 

We  are  performing  calculations  of  dissociative  attachment  of  low  energy  (<  5eV)  electrons 
to  molecular  hydrogen  in  vibrational  states  v  =  0  to  v  =  9,  using  a  nonlocal  complex 
potential.  Knowledge  of  these  cross  sections  is  important  for  optimizing  the  operating 
conditions  of  hydrogen  plasmas  from  which  H’  will  be  extracted.  The  present  work  is 
particularly  concerned  with  the  difference  between  the  rigorous  formulation  using  nonlocal 
potentials  (that  is,  integro-differential  equations)  and  the  simpler  approximation  assuming 
local  potentials  (differential  equations). 

A  manuscript  describing  the  calculations  has  been  accepted  for  publication  and  is 
attached  as  Appendix  E.  The  calculations  include  the  nonlocal  and  energy-dependent  width 
and  shift  operators  that  are  defined  by  resonant  scattering  theory. 
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We  compared  the  use  of  local  and  nonlocal  potentials  by  performing  calculations 
with  several  approximate  forms  of  the  full  nonlocal  potential.  We  found  that  in  certain 
cases  we  could  define  an  appropriate  local  potential,  but  only  by  first  calculating  and  then 
modifying  the  full  nonlocal  coupling  terms.  This  result  explains  the  relative  accuracy  of 
earlier  work  based  on  fitting  experimental  data  using  local  potentials,  because  it  shows  that 
one  does  expect  an  approximate  local  potential  to  exist.  Unfortunately,  the  appropriate 
local  potential  is  not  the  one  that  can  easily  be  identified  with  fixed-nuclei  properties  of  the 
resonant  electronic  state.  In  fact,  using  the  fixed-nuclei  curve  led  to  serious  errors.  Our 
conclusion  is  that  detailed  calculations  of  the  electronic  coupling  matrix  elements,  and 
inclusion  of  the  nonlocal  width  and  shift,  are  necessary  for  definitive  calculations  of 
electron  attachment. 
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APPENDIX  A 


Resonances  in  Dissociative  Recombination  of  Electrons  with  h£ 

From  Dissociative  Recombination:  Theory,  Experiment,  and  Applications 
Edited  by  J.B.A.  Mitchell  and  S.  Guberman  (World  Scientific  Publishing  Co.,  1988). 


Resonances  in  Dissociative  Recombination  of  Electrons  with  H2+ 


A.  P.  Hickman 

Molecular  Physics  Department 
SRI  International 
Menlo  Park,  CA  94025 


ABSTRACT 

Dissociative  recombination  of  electrons  with  diatomic 
molecular  ions  is  reviewed,  with  emphasis  on  the  qualita¬ 
tive  description  of  the  "direct"  and  "indirect"  mechanisms 
involved.  Calculations  for  H£+  are  presented  and  compared 
with  the  recent  experimental  results  of  Hus  et  al.  The 
calculations  confirm  the  importance  of  the  "indirect" 
mechanism,  which  accounts  for  the  resonant  structure 
observed  in  the  cross  section. 


1 .  INTRODUCTION 

Dissociative  recombination  (DR)  involves  the  capture  of  an 
electron  by  a  molecular  ion,  the  formation  of  a  compound  state,  and 
the  subsequent  dissociation  into  fragments.  This  paper  discusses  the 
mechanisms  that  come  into  play  as  the  process  unfolds  and  presents 
specific  calculations  for  e  +  H2+.  This  system  is  an  important 
prototype  because  it  involves  the  simplest  molecular  ion,  because 
experimental  data  ’  is  available,  and  because  the  relevant  excited 

a 

state  potential  curves  of  H2  have  been  calculated. 

Section  2  discusses  the  so-called  "direct"  and  "indirect" 
mechanisms.  Section  3  summarizes  the  theory.  A  single  resonance  in 
the  DR  cross  section  may  be  described  by  a  Beutler-Fano  profile,^  but 
the  general  approach  treats  overlapping  resonances.  Section  4  pre 
sents  calculations  for  H2+.  Section  5  contains  concluding  remarks. 
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2.  MECHANISMS  FOR  DISSOCIATIVE  RECOMBINATION 
2.1  Direct  Recombination 

The  physical  picture  of  direct  DR,  treated  by  Bardsley6  in  1968, 
may  be  understood  using  the  potential  curves  sketched  in  Fig.  la  and 
the  following  symbolic  representation: 

e  +  AB+  AB**  -►  A  +  B*  (1) 

After  capture  into  the  dissociating  state  AB  ,  the  nuclei  follow  the 
dissociating  potential  curve  V  .  The  excess  electronic  energy  of  the 
resonant  state  is  converted  to  kinetic  energy  of  the  dissociating 
products.  These  products  must  separate  to  the  point  where  V  crosses 
V+  for  the  recombination  to  be  complete.  Competing  with  the 
stabilization  of  the  resonant  state  by  dissociation  is  electronic 
autoionization  [indicated  by  the  leftward  pointing  arrow  in  Eq. 

(1)].  For  H2,  classical  mechanics  predicts  that  a  time  on  the  order 
of  10' 15  sec  is  required  for  dissociation.  Typical  calculated 
values^  of  the  autoionization  lifetime  are  also  on  the  order  of  10'^ 
sec.  The  conclusion  is  that  dissociation  provides  a  stabilization 
mechanism  fast  enough  to  compete  effectively  with  autoionization. 


RA-554SM 


Figure  1.  Illustration  of  direct  (a)  and  indirect  (b)  recombination. 
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The  details  of  the  theoretical  analysis  of  direct  DR  are  avail¬ 
able  in  the  literature . The  result  is  that  the  cross  section  is 
approximately  proportional  to  the  square  of  the  following  integral: 


“  (ir)1/2  f  *v_<R)  Vel(R)  F(R)  dR 


(2) 


'0  0 

is  the  wave  function  of  the  initial  vibrational  level  of  the 
molecular  ion,  is  the  matrix  element  for  electron  capture  (or 
autoionization) ,  and  F(R)  is  the  regular  solution  to  the  Schroedinger 
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equation  for  the  final  state  potential  V  (R)  at  energy  E  -  k  /2M. 

This  integral  is  essentially  a  Franck-Condon  factor  between  the  initial 
and  final  states,  mediated  by  the  electronic  coupling  V 

Direct  DR  for  thermal  energy  electrons  is  expected  to  have  a 
large  cross  section  when  the  Franck-Condon  factor  is  large.  This 
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situation  usually  occurs  when  the  dissociating  potential  V  crosses 
the  ionic  potential  V+  near  the  minimum  of  the  latter.  Conversely, 
when  the  dissociating  curve  does  not  cross  V+  near  the  minimum,  one 
expects  the  thermal  energy  cross  section  for  direct  DR  to  be  small. 
However,  even  when  direct  DR  is  small,  other  mechanisms  may  be  sig¬ 
nificant  and  may  lead  to  a  large  total  cross  section.  A  qualitative 
discussion  of  this  possibility  is  presented  in  the  next  subsection. 


2.2  Indirect  Recombination 


At  certain  well-defined  energies,  the  incident  electron  may  be 
captured  into  compound  resonant  states,  here  denoted  AB  (n,v).  These 

•ft  »»*•_ 

states  provide  an  alternative  pathway  to  the  final  state  A  +  B  .  The 
situation  represented  by  Eq.(l)  must  therefore  be  generalized: 


e  +  AB 


- AB 

#  X 

AB*(n,vr 


->►  A  +  B 


(3) 


The  additional  states  are  labelled  by  the  indices  n  and  v.  Their 
physical  form  is  an  ion  core  AB+  in  vibrational  level  v,  and  a 
Rydberg  electron  in  an  orbital  whose  principle  quantum  number  is  n. 
(We  omit  discussion  of  the  electron's  angular  momentum). 
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One  way  to  consider  these  resonant  states  is  the  following. 
Beginning  with  the  ionic  potential  V+(R) ,  define  the  series  of 
Rydberg  potential  energy  curves  Vn(R)  -  V+(R)  -  0.5  [n-/i(R)] The 
quantum  defect  /i(R)  defines  the  binding  energy  of  the  Rydberg  elec¬ 
tron.  Each  potential  Vn(R)  of  the  molecule  AB  has  a  set  of  vibra¬ 
tional  energy  levels.  One  of  the  Vn(R)  is  shown  in  Fig.  lb,  and  one 
of  its  vibrational  energy  levels  is  drawn.  For  an  incident  electron 
of  the  proper  energy,  the  total  energy  of  the  system  will  be  close  to 
that  of  this  particular  AB  (n,v) .  One  expects  the  cross  section  for 
recombination  to  have  additional  structure  near  this  energy. 

3 .  THEORY 

6  8 

The  theory  of  DR  may  be  formulated  using  the  matrix  elements 

Vv  “  r  S  V(R,)  Vel(R,)  G(R'’R>  Vel(R)  *V(R)  dR' 

where  G(R' ,R) ,  defined  explicitly  in  ref.  8,  is  a  Green  function  for 
the  Schroedinger  equation  with  potential  V  .  The  matrix  element 
is  related  to  the  amplitude  for  capture  into  the  resonant  state 
AB#(n,v')  from  the  initial  vibrational  state  v  of  AB+.  The  mechanism 
for  this  capture  is  a  second  order  (in  Vg^)  electronic  coupling  via 
the  resonant  state  AB  .  The  momentum  transfer  from  the  light 
electron  to  the  heavy  nuclei  via  potential  scattering  is  neglected. 
This  approximation  is  well-fulfilled  in  practice  for  H2+. 

For  the  present  discussion,  it  is  sufficient  to  note  that  G 
corresponds  to  outgoing-wave  boundary  conditions  (i.e.,  dissociation) 
and  is  therefore  a  complex-valued  function.  We  may  write 

gv,v  -  Re  gv,v  +  <i/*Kv»£v  (4) 

where  the  £v  are  defined  in  Eq.  (2). 

Giusti^  first  treated  DR  including  both  direct  and  indirect 
mechanisms,  using  "first  order"  Multichannel  Quantum  Defect  Theory 
(MQDT) .  The  terms  Re  gy«v  were  neglected.  Later  calculations  for 
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e  +  i*2+  by  Giusti  et  al.^  and  Nakashima  et  al.^®  also  neglected 
Re  g v»v.  Hickman  included  Re  gy>v  in  an  approach  based  on  the 
Independent  Resonance  Approximation  (IRA)  and  discovered  that  the 
previously  neglected  terms  were  quite  important  for  e  +  H2+. 
However,  the  relation  between  the  two  theories  (MQDT  and  IRA)  was 
obscure.  Our  recent  work,  summarized  here,  clarifies  this  relation 
and  ;  .fies  the  two  theories. 

MQDT  is  formulated  in  terms  of  a  reactance  matrix  coupling 
several  different  arrangement  channels .  These  channels  have  the 

i 

asymptotic  forms  A  +  B  (labelled  d) ,  and  AB  (v)  +  e,  for  v-0,1,... 
The  specific  form  of  the  matrix  that  we  use  is 

d  v-0  v-1  v-2  v-3 


(5) 


This  matrix  may  be  "contracted"  into  an  S  matrix  using  the  general 
MQDT  prescription.^  S  may  have  a  smaller  dimension  than  the  reac¬ 
tance  matrix,  because  some  of  the  arrangement  channels  are  closed. 

We  have  found  a  limiting  case  in  which  the  MQDT  contraction  can 
be  explicitly  compared  with  the  IRA.  Consider  the  channels  d,  v-0 
(the  initial  channel),  and  v-1  (a  closed  channel).  The  effect  of  a 
single  resonance  on  the  recombination  may  then  be  described  in  terms 
of  a  background  process  (direct  DR)  that  is  modified  in  a  restricted 
energy  range  by  the  resonance.  Identical  result  are  obtained  with 
the  IRA  and  MQDT: 

.  .  direct.  .  .  .  .2,,  2  , 

<7(0  -  a  ( O  (x  +  q)  /(x  +  1)  (6) 

where  x  -  (E-Eres)/(r/2) ,  and  T  -  2£^/jm^  is  the  width  of  the 
resonance.  We  find  that  q  is  given  by  -(Re  gQ^)/(Im  ggi)*  This 
result  demonstrates  the  equivalence  of  MQDT  and  the  IRA  for  single 
resonances ,  provided  that  one  includes  the  effect  of  the  full  Green 
function  as  shown  in  Eq.  (5).  The  "first-order"  limit  of  MQDT 

Q 

discussed  by  Giusti  corresponds  to  neglecting  the  block  of  terms  in 
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Eq.  (5)  whose  values  are  w  Re  gy*v.  For  the  calculations  reported 
here,  we  use  the  full  MQDT  approach  based  on  the  reactance  matrix, 

Eq.  (5). 

4.  RESULTS  FOR  e  +  H2+ 

Fig.  2  shows  the  intermediate  states  AB*(n,v)  introduced  in 
2.2.  The  vibrational  levels  of  the  ion  AB+  are  drawn,  as  well  as  an 
electronic  Rydberg  series  converging  to  each  ionization  limit. 
Clearly,  for  many  incident  electron  energies,  contributions  from  more 
than  one  intermediate  state  must  be  included. 

Two  limiting  cases  may  be  identified:  small  Av  and  large  n  (the 
left  hand  side  of  Fig.  2),  and  large  Av  and  small  n  (the  right  hand 
side  of  Fig.  2).  The  discussion  in  2.2  and  Fig.  lb  correspond  to 
resonances  of  the  first  type.  One  may  wonder  whether  resonances  of 
the  second  type,  involving  up  to  seven  quanta  of  vibrational  energy 
transfer,  are  important.  Following  the  arguments  of  O'Malley,  one 


RA-5549-2 

Figure  2.  Rydberg  energy  levels  of  H2.  Inset:  Curve  crossing 
mechanism  for  high  v,  low  n  resonances. 
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sees  that  such  resonances  may  be  interpreted  in  terms  of  a  curve 
crossing  as  shown  in  the  upper  left  of  Fig.  2.  The  electron  is  cap¬ 
tured  into  the  dissociating  state,  but  during  dissociation  the  system 
undergoes  a  curve  crossing  to  the  potential  curve  Vn(R) ,  where  n-3 
(the  smallest  value  of  n  possible  for  the  da  waves  considered  here). 

The  effect  of  the  overlapping  resonances  may  be  probed  by  model 
calculations  that  include  only  selected  resonances.  Two  series  of 
such  calculations  are  presented  in  Fig.  3.  First,  we  consider  the 
case  Vq-0,  (vq  denotes  the  initial  vibrational  level  of  H2+) .  In 
this  case  the  Franck- Condon  factor  is  low  and  one  expects  the  direct 
mechanism  to  be  small.  Panel  (a)  shows  a  dotted  line  corresponding 
to  the  direct  process;  the  transition  probability  is  indeed  small, 
about  10%.  This  panel  also  shows  the  regular  structure  arising  from 
the  Rydberg  series  converging  to  v-1  (cf.  Fig.  2).  The  structure 
defined  by  Eq.  (7),  which  could  be  obtained  by  the  IRA,  is  essen¬ 
tially  repeated,  properly  scaled,  by  each  member  of  the  Rydberg 
series.  Panel  (b)  shows  the  result  if  one  considers  a  single 
resonance  corresponding  to  v«7  and  n-3.  Panel  (c)  shows  the  combined 
results.  The  overall  pattern  is  more  complicated,  but  still 
reminiscent  of  the  simpler  results  obtained  in  (a)  and  (b) .  Note 
that  the  deep  minimum  seen  in  (b)  near  0.20  eV  persists. 

Panels  (d)-(f)  show  a  parallel  calculation  for  the  case  vq-1. 

In  this  case  the  Franck-Condon  factor  is  large,  and  the  direct 
mechanism  (dotted  line)  is  much  larger,  about  70%.  The  series  of 
resonances  for  v-1  to  v-2  excitation  [panel  (d)]  are  "window" 
resonances,  that  is,  purely  destructive  interference.  The  single 
resonance  for  v-1  to  v-8  excitation  is  also  of  this  type.  Again,  the 
composite  case  [panel  (f)]  exhibits  features  that  are  similar  to 
those  of  each  contributing  resonance. 

Our  results  suggest  that  the  indirect  mechanism  tends  to  inhibit 
recombination  when  the  direct  mechanism  is  large.  When  the  direct 
mechanism  is  small,  however,  the  indirect  mechanism  may  enhance  or 
inhibit  the  cross  section.  Having  a  smaller  cross  section  when  the 
indirect  mechanism  is  included  is  consistent  with  the  physical  picture 
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developed  by  O'Malley  based  on  classical  trajectories.  The  curve 
crossings  in  the  outgoing  channel  make  dissociation  more  difficult. 

The  exact  energies  of  the  resonances  corresponding  to  Urge  Av 
are  sensitive  to  the  molecular  potential  curves.  We  find  thdt  dif¬ 
ferent  values  of  /i(R)  from  the  literature  lead  to  different  values  of 
the  resonance  energies.  For  Av  -  7-8,  the  discrepancies  are  large, 
on  the  order  of  tenths  of  eV.  (Another  contribution  to  the  resonance 
energy  comes  from  Re  g^) .  Since  the  effect  of  these  large  v 
resonances  is  often  to  reduce  substantially  the  cross  section  over  a 
particular  energy  range,  it  is  essential  to  calculate  accurately  the 
location  of  each  resonance. 

These  uncertainties  preclude  a  definitive  calculation  even 
though  the  MQDT  formalism  provides  the  computational  framework. 
Nevertheless,  Fig.  4  presents  calculations  including  the  vibrational 
levels  v-0-9  for  H2+.  Separate  calculations  for  Vq-0  and  Vq-1  were 
performed;  the  experimental  values  of  Hus  et  al.,  shown  for 
comparison  in  both  panels,  correspond  to  a  mixture  of  these  two 
initial  states.  Some  testing  of  the  effects  of  using  different 
functions  /*(R)  from  the  literature  was  carried  out.  Our  results  tend 
to  favor  the  /i( R)  of  Nakashima  et  al.,*®  which  is  numerically  much 
smaller  than  the  other  functions.  Those  results  are  shown. 
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Figure  4.  Comparison  with  experiment  (Hus  et  al.  )  for  e  +  H2. 
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5 .  CONCLUDING  REMARKS 

There  has  been  considerable  progress  in  the  understanding  of  DR 
of  molecular  ions  in  the  past  few  years.  The  recent  experimental 
studies  of  Hus  et  al.  have  an  electron  energy  resolution  of  about 
0.02  eV,  and  reveal  a  complex  structure  of  the  cross  section. 
Theoretical  calculations  confirm  that  the  structure  arises  from  the 
pattern  of  overlapping  Rydberg  series  converging  to  the  various 
vibrational  levels  of  the  target  molecular  ion.  These  Rydberg  states 
collectively  contribute  to  the  "indirect"  mechanism  of  recombination, 
which  is  of  prime  importance  for  e  +  H2+.  This  indirect  mechanism 
should  be  assessed  in  future  calculations  on  other  systems, 
particularly  in  cases  where  Franck- Condon  arguments  suggest  that  the 
direct  mechanism  will  be  small. 
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ABSTRACT. 

High  resolution  measurements  for  the  dissociative 
recombination  of  H2*  (v=0)  are  presented  and  compared  with 
theoretical  calculations.  The  observed  structure  is  explained  with 
a  model  that  invokes  the  coupling  to  intermediate  electronic  states 
belonging  to  two  separate  Rydberg  series. 


INTRODUCTION. 

The  complex  collision  process  of  dissociative  recombination 
(DR)  can  be  better  understood  by  comparing  high  quality 
experimental  results  with  theoretical  models.  The  most  likely 
candidate  for  such  a  comparison  is  H2*,  the  simplest  molecular  ion. 
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For  this  system,  DR  can  be  schematically  written 


e  +  H2*  -»  H  +  H*  ...  (4) 

A  number  of  theoretical  calculations  have  been  performed  for  this 
ion1"4.  These  calculations  have  predicted  that  the  cross  section 
for  electrons  recombining  with  H2+  in  a  specific  vibrational  state 
should  display  a  series  of  narrow  resonances  due  to  the  influence 
of  the  neutral  Rydberg  states  lying  below  the  ion  state.  In  order 
to  observe  these  resonances  it  is  necessary  to  prepare  the  ions  in 
a  single  vibrational  state  and  to  have  a  collision  at  a  well 
defined  energy.  The  most  suitable  method  for  doing  this  is  to  use 
the  merged  beam  technique  which  has  been  developed  at  the 
University  of  Western  Ontario.  A  preliminary  measurement  that 
demonstrated  the  existence  of  the  resonances  has  already  been 
published5.  This  paper  describes  some  recent  measurements  that 
were  performed  using  very  high  energy  resolution  and  a  fine  energy 
mesh.  The  results  are  now  of  sufficient  accuracy  to  associate 
particular  features  of  the  cross  section  with  specific  Rydberg 
states.  Additional  theoretical  work  is  necessary  to  compare 
experiment  and  theory  at  this  level  of  accuracy.  Extensions  to 
previous  theory4,  and  the  very  good  comparison  with  experiment 
thereby  obtained,  are  also  presented  in  this  paper. 

EXPERIMENTAL  METHOD. 

The  Merged  Beam  Electron  Ion  Beam  Experiment  at  the  University 
of  Western  Ontario  has  been  described  in  detail  elsewhere6  and  only 
a  brief  overview  of  the  apparatus  will  be  presented  here.  A  beam  of 
H2+  ions  are  produced  by  electron  impact  in  a  radiofrequency  ion 
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trap  source  located  in  the  terminal  of  a  Van  de  Graaff  accelerator. 
This  source,  when  operated  with  a  low  extraction  potential  is 
capable  of  producing  ions  with  very  low  internal  energies5,7.  The 
ions  are  accelerated  to  400  kev,  mass  analyzed  and  passed  into  the 
interaction  region  where  they  are  merged  with  a  beam  of  electrons 
erived  from  an  indirectly-heated,  barium  oxide  cathode,  using  a 
trochoidal  analyzer.  The  two  beams  interact  over  a  distance  of  8.6 
cm  before  being  demerged.  The  electron  beam  is  collected  in  a 
Faraday  cup,  the  ion  beam  is  deflected  electrostatically  into  a 
second  Faraday  cup  and  the  undeflected  neutrals  impinge  upon  a 
surface  barrier  detector  and  are  detected.  The  electron  beam  is 
modulated  and  the  signals  are  measured  in  and  out  of  phase  with  the 
modulation  in  order  to  distinguish  true  recombination  signals  from 
neutrals  produced  in  collisions  between  the  ion  beam  and  the 
background  gas.  The  overlap  of  the  two  beams  is  measured  at  two 
places  in  the  interaction  region  and  the  effective  collision  area, 
F,  is  determined. 

The  cross  section  for  recombination  is  determined  using  the 
formula 


o- 


cnea  I 

IlI*L  1^-vJ 


.  (1) 


where  Cn  is  the  neutral  count  rate,  e,  the  electronic  charge,  I4, 
v, ,  Ie,  and  ve  are  the  ion  and  electron  beam  currents  and  velocities 
respectively,  L  is  the  length  of  the  interaction  region  and  F  is 
the  effective  collision  area. 

The  energy  in  the  center-of-mass  is  given  by: 

Ec*  =  (EeVZ  "  E+1/2)2  +  Ee1/2Ej/2  02  ...(2) 
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where  8  is  the  intersection  angle  between  the  two  beams,  Ee  is  the 
electron  beam  energy  and  E*  the  reduced  ion  energy  is  given  by  E4 
=  Ej.me/mj,  Ej  being  the  ion  beam  energy.  By  differentiating  this 
equation  with  respect  to  E+,  Ee  and  0,  it  is  possible  to  write  down 
an  expression  for  the  uncertainty  in  the  center-of-mass  in  this 
measurement . 

AE^-K [  (1-  ( )  1/23  A Ej*+< [1-  (^  > 1/2]  aeJ2+ 

+  [2(E8EJ1/20A0]2]1/2.  .  .  (3) 

It  can  be  seen  from  equation  (3)  that  when  Ee  =  E+,  the  first 
two  terms  go  to  zero  and  the  uncert* i:.cy  is  determined  mainly  by 
the  uncertainty  in  the  intersection  angle.  Using  well  collimated 
beams  and  careful  alignment  this  can  be  made  quite  small.  The 
unique  property  of  the  merged  beam  technique  is  that  the  ultimate 
energy  resolution  at  low  energies  is  rather  insensitive  to  spreads 
in  the  energies  of  the  electron  and  ior  beams. 

The  energy  scale  of  the  measurement  was  calibrated  using  a 
method  described  in  ref.8  in  which  the  ion  beam  energy  is  first 
determined  using  a  nuclear  resonance  technique  and  the  electron 
energy  is  subsequently  obtained  from  the  position  of  the  maximum  in 
the  DR  signal  as  the  electron  energy  is  scanned  through  the  Ee«E+ 
energy  region.  E^  can  then  be  calculated  using  Eq.2.  The  largest 
uncertainty  associated  with  this  calibration  is  in  the 
determination  of  E#  and  this  amounts  to  ±  0.5  eV.  From  Eq.2  it  can 
be  seen  that  this  inaccuracy  leads  to  an  uncertainty  in  the  center- 
of-mass  energy  of  ±  7meV  at  20  meV  andil5  meV  at  100  meV. 
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RESULTS  AND  DISCUSSION. 


The  cross  sections  for  the  dissociative  recombination  of  H2* 
ions  with  low  internal  energy,  measured  in  the  present  study  are 
shown  in  figure  la.  The  ion  and  electron  beam  currents  used  were 
20jiA  and  lxlO'uA  respectively  and  more  than  24  hours  per  point  were 
needed  in  order  to  achieve  the  level  of  statistical  accuracy  shown 
here.  The  error  bars  shown  represent  one  standard  deviation  and 
systematic  errors,  mainly  associated  with  measurement  of  F, provide 
an  uncertainty  in  the  absolute  cross  sections  of  ±  20%.  The  effect 
of  systematic  errors  on  the  relative  uncertainty  in  the  data  points 
is  about  4%  however  since  F  is  a  quantity  which  changes  slowly  with 
energy. 

The  merged  beam  technique  is  very  stable  regarding  the 
collision  energy  and  measured  points  are  very  repeatable.  Three 
well  defined  window  resonances  are  clearly  apparent  in  the  data. 
From  the  figure  it  can  be  seen  that  at  the  low  end  of  the  energy 
scale,  the  resolution  must  be  of  the  order  of  5  mev.  Using 
equation  (3)  this  translates  to  a  value  of  about  10  meV  at  the 
upper  end.  At  a  higher  extraction  potential  the  central  resonance 
was  found  to  disappear  and  this  is  due  the  population  of  v=l  states 
in  the  beam  increasing  causing  this  resonance  to  be  washed  out.  By 
comparison  with  theoretical  results  it  is  believed  that  the  results 
presented  in  figure  la  refer  only  to  v=0  ions.  This  is  discussed 
further  below.  No  effort  was  made  to  restrict  the  rotational 
population  of  the  beam  and  the  source  gas  was  undoubtedly  a  mixture 
of  ortho  and  para  hydrogen.  The  rotational  population  must 
therefore  be  reflect  a  thermal  distribution  of  odd  and  even  J 
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levels. 


The  direct  and  indirect  mechanisms  of  dissociative 
recombination  have  been  identified  and  discussed  in  the  literature 
1-4 .  One  can  augment  Eq.  (1)  as  follows 

e  +  H2*  -  H2**  -  H  +  H*  ...(4) 

(v,n, 1,A) 

The  top  line  of  Eq. (4)  corresponds  to  direct  recombination:  an 
incident  electron  is  captured  directly  into  the  dissociating 
state  and  the  nuclear  motion  on  this  potential  leads  to  the  final 
products.  At  certain  well-defined  energies,  the  electron  may  also 
be  captured  into  compound  resonance  states,  here  denoted  on  the 
second  line  by  H2(v,n,l,A).  The  physical  form  of  such  a  state  is 
an  ion  core  H2+  in  vibrational  level  v  and  a  Rydberg  electron  whose 
quantum  numbers  are  n,  1  and  A  (projection  on  the  internuclear 
axis)  .  This  state  is  formed  when  an  electron,  incident  on  H24  (v=0) 
vibrationally  excites  the  ion  to  v=l,  thereby  losing  energy  and 
dropping  down  into  a  Rydberg  orbital.  The  state  H2(v,n,l,A)  may 
decay  to  H  +  H*  by  predissociation.  This  sequence  of  events 
provides  an  alternative  mechanism  for  recombination,  and  leads  to 
classic  interference  phenomena.  The  direct  mechanism  provides  a 
slowly  varying  background  cross  section;  at  those  energies  for 
which  the  intermediate  compound  state  is  allowed,  the  cross  section 
exhibits  rapidly  varying  resonance  structure.  By  analyzing  this 
structure,  we  can  infer  information  about  the  compound  states. 

The  specific  resonant  states  H2(v,n,l,A  involved  have  1=0  or 
2  and  A=  0,  since  the  coupling  is  through  the  's  4  (predominantly 
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2p a2)  state.  These  resonances  correspond  to  sa  and  do  Rydberg 
series.  Our  analysis  indicates  that  both  series  must  be  included. 
We  do  this  by  generalizing  the  theoretical  approach  used  in  ref.  4. 
The  matrix  element  leading  to  direct  recombination  was  written 

lv-(lf-)1/2{xv(R)  -  Vb1{R)  .F(R)  .dR...  (5) 

where  M  is  the  reduced  mass  of  H2,  xv(R)  is  wave  function  of  the 
vth  vibrational  level  of  H2*,  Vel(R)  is  the  matrix  element  for 
electron  capture  (or  autoionization),  F(R)  is  the  regular  solution 
to  the  Schrodinger  equation  for  the  final  state  potential  V*(R)  at 
energy  E=*i2k2/2M .  This  integral  is  essentially  a  Franck-Condon 
factor  between  initial  and  final  states,  mediated  by  the  electronic 
coupling  Vel . 

The  indirect  mechanism  is  introduced  using  the  matrix  element 
9v'vm{Xv'(R' >  >Vel(R')  . GiR',R )  ,V0l(R)  .XV(R)  .  dR'dR.  .  .  (6) 


where  G(R',R) ,  defined  explicitly  in  ref. 3,  is  a  Green  function  for 
the  Schrddinger  equation  with  potential  V*(R).  The  matrix  element 
gv,v  is  related  to  the  amplitude  for  vibrational  excitation  of  H2* 
from  state  v  to  v'  by  an  incident  electron.  The  two  vibrational 
states  are  not  directly  coupled  to  each  other,  but  each  is  coupled 
electronically  to  the  resonant  state  H2**. 

The  coupling  of  the  resonant  state  H2**  to  both  sa  and  da 
electrons  is  treated  in  an  approximate  way  by  assuming  that  the 
ratio  of  the  couplings  to  the  two  partial  waves  is  a  constant. 
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independent  of  R.  Then  we  can  extract  the  necessary  matrix 
elements  from  available  literature  values  of  the  total  widths  by 
replacing  the  coupling  matrix  element  Vel(R)  by  c^V^fR)  (1=0  for 
coupling  to  so  and  1=2  for  do).  Then  a02  +a22  =1,  and  we  can  write 

•  ■  (7) 


’  1  1 

We  use  a  theoretical  approach  based  on  Multichannel  Quantum 
Defect  Theory  (MQDT) ,  which  was  formulated  by  Giusti9  and  modified 
by  Hickman4  to  include  higher  order  effects.  MQDT  is  formulated  in 
terms  of  a  reactance  matrix  K  coupling  to  different  arrangement 
channels.  These  channels  have  the  asymptotic  form  H  +  H*  (labelled 
d)  and  H2*(v)  +  e(lo)  (v=0,l,...  and  1=0,2).  The  rotational  state 
of  the  H2*  ion  has  not  been  included.  The  specific  states  we 
include  in  the  present  calculation  are  therefore  (d) ,  (so,  v=0)  , 
(do,  v=0) ,  (so,  v=l) ,  and  (do,  v=l) ,  and  the  K  matrix  is 


K  = 
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...(9) 


The  diagonal  blocks  contain  the  phase  shifts  and  nd,  which  are 
related  to  the  energies  of  the  autoionizing  states  H2(v=l,n, 1,A) . 
This  matrix  is  "contracted  "  into  an  S  matrix  for  DR  using  the 
general  MQDT  prescription10.  S  has  a  smaller  dimension  that  the  K 
matrix  for  the  energies  of  interest,  for  which  the  v=l  channels  are 
closed. 

The  K  matrix  given  in  equation  (6)  provides  flexibility  to 
model  the  data.  The  quantum  defects  and  nd  determine  the 
positions  of  the  two  series  of  resonances,  given  by 


e-A£- - ± - 

2  (fl-Hj)2 


(8) 


where  AE  is  the  vibrational  excitation  energy  from  v=0  to  v=l,  1=0 
and  2  and  n  takes  on  integral  values.  By  inverting  equation  (7) , 
we  can  estimate  the  effective  values  n*=n-Mt  of  the  resonances 
observed  experimentally.  The  three  resonances  shown  in  figure  1 
correspond  to  n*=7 . 3 ,  7.7  and  8.1.  Taking  into  account  the 
uncertainty  in  the  energy  scale  of  the  data,  we  can  approximately 
model  the  data  by  assuming  that  the  middle  resonance  (n*  =7.7) 
belongs  to  one  series  and  the  other  two  (n*  =7.3,  8.1)  to  the  other 
series. 

Calculations  have  been  performed  to  investigate  these 
alternatives.  The  values  of  the  parameters  in  equation  6  were 
estimated  from  the  literature  and  then  systematically  varied  to 
achieve  a  good  fit.  Hazi  et  al.11  determined  a0  =  0.86  and  a2  =  0.54. 
The  quantum  defects  can  be  taken  from  ref.11  or  directly  estimated 
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to  be  jig=-0 . 108  and  nd  =0.053  from  the  H2  GK  or  HHfbarJ  potential 
curves.  These  values  can  only  be  considered  guidelines;  adjusting 
them  is  not  unreasonable.  The  quantum  defects  are  expected  to  vary 
weakly  with  energy,  and  their  value  the  ^ 4  continuum  state  will 
certainly  differ  from  their  value  in  the  region  of  the  GK  and 
HH[bar]  states. 

Determining  which  resonances  are  sa  and  which  are  da  is  not 
possible  on  the  basis  of  the  quality  of  the  fit  to  the  data.  Both 
alternatives  lead  to  a  plausible  fit  for  the  appropriate  values  of 
the  quantum  defects  /us  and  nd  and  the  branching  ratios  a0  and  a2. 
)n  physical  grounds,  one  generally  expects  the  d  quantum  defect  to 
be  smaller  than  the  s.  Also  Hazi  et  al.11  considered  a  fit  to  the 
EF  and  GK  potential  curves  of  H2  and  concluded  that  the  width  of 
the  sa  partial  wave  is  smaller  than  the  da.  Based  on  these 
considerations,  we  favour  the  assignment  that  the  45  meV  resonance 
is  an  sa  state  and  the  other  two  are  da  states.  For  this 
assignment,  we  used  ms=0.4  and  ^d=-0.1. 

Figure  lb  shows  the  comparison  calculated  DR  cross  section. 
The  model  based  on  two  rydberg  series  provides  an  excellent 
interpretation  of  the  data.  Furthermore,  the  data  cannot  be  fit  by 
assuming  a  single  rydberg  series. 
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FIGURE  CAPTIONS. 


1.  Experimental  (a)  and  theoretical  (b)  cross  sections  for  the 
dissociative  recombination  of  H2*(v=0). 
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Two-photon  excitation  and  excited-state  absorption  cross  sections 
for  H2  E,F  (v  =6):  Measurement  and  calculations 
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Douglas  J.  Bamford/  and  William  K.  Bischel 
Molecular  Physics  Laboratory,  SRI  International,  Menlo  Park,  California  94025 
(Received  21  November  1988) 

The  absolute  two-photon  excitation  cross  section  for  the  H2  E,F  '2,  (u’=6)«— Y  '2,  (u"=0)  Q(l) 
transition  at  193  nm  has  been  measured  by  observing  the  E,F  '2*  ( t> '  =  6 ) — '24  (u"=0)  fluores¬ 
cence  at  ~750nm.  The  measured  integrated  two-photon  excitation  cross  section,  (2.0±0.9)X  10“ ,6 
cm4,  is  in  good  agreement  with  the  theoretical  value  of  2.8 X  10“ 36  cm4,  which  is  obtained  from  pre¬ 
viously  published  calculations  [Huo  and  Jaffe,  Chem.  Phys.  Lett.  101,  463  (1983)].  The  absolute 
cross  section  for  photoabsorption  by  the  E,F  state  at  355  and  193  nm  was  also  measured  by  moni¬ 
toring  the  depletion  of  the  750-nm  fluorescence  caused  by  a  second  laser.  The  measured  cross  sec¬ 
tions  for  photoabsorption  by  the  E,F  (u=6)  state  are  (9.7±2.4)X  10'18  cm2  at  355  nm  and 
(6.4±  1 . 3 ) X  I0-18  cm2  at  193  nm.  Comparison  with  theoretical  estimates  [Cohn,  J.  Chem.  Phys,  57, 
2456  (1972)]  cf  the  direct  photoionization  cross  section  of  this  state  (8.5X  1CT'8  cm2  at  355  nm  and 
3.2 X  10“ 18  cm2  at  193  nm)  suggests  that  other  processes  may  contribute  to  the  photoabsorption. 
Photoexcitation  of  the  dissociating  '2„  autoionizing  state  is  found  to  be  important.  Cross  sections 
for  this  additional  channel,  which  may  lead  to  dissociation  (H*  +H  or  H*  +H~)  or  auloionization 
(H,*  +e ” ),  are  calculated  to  be  3.6X  10~'8  cm2  at  355  nm  and  7.6X  10" 18  cm2  at  193  nm.  The  cal¬ 
culated  branching  ratio  strongly  favors  dissociation  at  355  nm  and  autoionization  at  193  nm. 


I.  INTRODUCTION 

Molecular  hydrogen  occupies  a  position  of  unique  im¬ 
portance  in  chemical  physics  because  it  is  the  simplest 
neutral  molecule.  As  such,  it  supports  a  rich  interplay 
between  experiment  and  theory,  due  largely  to  the  rela¬ 
tive  tractability  of  large-scale  calculations  on  this  four- 
particle  molecular  system.  Hydrogen  is  also  of  consider¬ 
able  practical  importance  in  areas  such  as  combustion, 1 
extraterrestrial  atmospheres,2  vacuum-uv  laser  sources,3 
laser-based  analytical  probes,4  and  numerous  others.  In 
view  of  this,  we  have  undertaken  a  series  of  experiments 
to  compare  absolute  cross-section  measurements  to  calcu¬ 
lations  for  two-photon  resonant  multiphoton  transitions 
in  Hj.  Ab  initio  calculations  of  the  two-photon  transition 
moment5  for  H2  E,F  (v'  =  (>)— X  (v"=0)  have  been  avail¬ 
able  in  the  literature  for  several  years  and  can  be  used  to 
compare  to  quantitative  measurements.  In  addition,  this 
research  was  motivated  by  the  need  to  develop  detection 
methods  for  vibrationally  excited  ground  electronic  state 
hydrogen  for  application  to  the  diagnostics  of  H_  plasma 
sources.6,7  We  expect  the  results  reported  here  for  the 
i)"  =  0  ground  vibrational  state  will  serve  as  a  benchmark 
for  determining  the  sensitivity  for  H2  (u”>0)  detection 
using  two-photon  excitation. 

In  this  paper  we  report  absolute  measurements  of  both 
the  two-photon  cross  section  for  excitation  in  the  transi¬ 
tion  from  the  X  (u"  =  0)  to  the  E,F  (v'  =  6)  state  of  H2, 
and  cross  sections  for  absorption  by  the  excited  state. 
The  wavelength  dependence  of  the  absorption  cross  sec¬ 
tion  is  of  particular  interest  since  several  mecha¬ 
nisms  in  addition  to  direct  photoionization  can  contrib- 


FIG.  I.  Selected  energy  levels  for  Hj.  The  threshold  for  ion¬ 
ization  is  124417  cm' for  dissociative  ionization  the  threshold 
is  145  796  cm  1  (Ref.  35). 
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ute.®-10  As  a  first  experiment  in  this  direction,  we  have 
measured  the  total  absorption  cross  section  at  two  wave¬ 
lengths  (193  and  355  nm).  These  results  are  compared  to 
theoretical  calculations  of  the  cross  sections  for  direct 
photoionization  and  for  other  channels  contributing  to 
absorption. 

The  H2  energy  levels  relevant  to  the  work  reported 
here  are  given  in  Fig.  1.  The  J'=  1  rotational  level  of  the 
E,F  (v'  =  6)  level  was  excited  using  two-photon  absorp¬ 
tion  at  193  nm.  The  excited  level  either  fluoresced  in  the 
near-ir  region  (750-840  nm)  or  absorbed  a  third  photon 
(either  at  193  or  355  nm).  The  experiment  consisted  of 
quantitatively  measuring  either  the  fluorescence  pro¬ 
duced  in  the  two-photon  process  or  the  depletion  of 
fluorescence  caused  by  absorption  of  a  third  photon,  and 
measuring  the  laser  intensity  in  the  observation  volume. 
Using  techniques  previously  developed  in  this  laborato¬ 
ry, 1112  two-photon  excitation  cross  sections  were  derived 
from  the  data.  Using  the  fluorescence  depletion  tech¬ 
nique,  15  excited-state  absorption  cross  sections  were  also 
obtained. 

The  paper  is  organized  as  follows.  Section  II  describes 
the  experiment,  including  a  detailed  description  of  two 
methods  used  to  produce  tunable  radiation  at  193  nm. 
Section  III  describes  the  theoretical  calculations.  Two- 
photon  excitation  cross  sections  for  the  E,F—X  transi¬ 
tion,  and  for  the  direct  photoionization  of  the  E,F  state, 
are  derived  from  work  already  in  the  literature.  Photoex¬ 
citation  to  a  dissociating,  autoionizing  state  from  the 
E,F  state  is  treated.  Section  IV  compares  and  discusses 
the  experimental  and  theoretical  results,  and  Sec.  V 
presents  concluding  remarks. 

II.  EXPERIMENT 
A.  Apparatus  and  procedure 

The  apparatus,  methods  of  data  acquisition,  and 
analysis  used  here  for  the  two-photon  absorption  cross- 
section  measurements  have  been  previously  described  in 
detail  for  atomic  oxygen.12  The  apparatus  is  diagramed 
in  Fig.  2.  Tunable  laser  radiation  at  193  nm  was  generat¬ 
ed  using  two  different  methods.  This  allowed  the  com¬ 
parison  of  cross-section  measurements  using  coherent 
sources  that  had  significantly  different  photon  statistics, 
bandwidths,  and  spatial  temporal  characteristics.  In  the 
first  method  a  frequency  doubled  Nd:yttrium-aluminum- 
gamet  (YAG)-pumped  dye  laser  (Quanta-Ray 
DCR2  +  PDL)  was  Raman  shifted  in  H2  to  generate  over 
100  j/J  ir.  the  fourth  anti-Stokes  (AS)  order  at  193  nm. 
Because  the  first  Stokes  order  of  the  Raman-shifted  dye 
laser  fundamental  has  nearly  the  same  wavelength  as  the 
fluorescence  being  observed  ( ~750  nm),  it  was  necessary 
to  prevent  any  stray  light  from  the  Raman  cell  from 
reaching  the  photomultiplier.  This  was  achieved  by 
focusing  the  193-nm  beam  (fourth  AS  order)  with  a  20- 
cm  lens  through  a  small  pinhole  in  an  opaque  enclosure 
which  surrounded  the  Raman  cell  output  and  the  Pellin- 
Broca  prism  used  to  separate  the  various  beams.  A 
second  20-cm  lens  outside  the  enclosure  then  reconverged 
the  193-nm  beam  to  a  focus  — 100  cm  away  in  the  experi- 


F1G.  2.  Apparatus  for  measurement  of  the  excited-state  ab¬ 
sorption  cross  section  at  193  nm  in  H;.  The  E,F  (u'  =  6)— X 
(i>"  =  0)  pfl)  two-photon  absorption  cross  section  was  measured 
in  a  similar  apparatus,  with  the  ArF  laser  beam  blocked,  and 
the  last  20-cm  lens  removed.  To  measure  the  excited-state 
single-photon  absorption  cross  section  at  355  nm,  the  potassium 
dihydrogen  phosphate  (K.DP)  crystal  was  used  to  generate  both 
532-nm  radiation  do  pump  the  dye  laser)  and  355-nm  radiation. 
The  latter  follows  a  path  like  that  of  the  ArF  laser  beam. 


mental  cell.  Pulse  energies  of  up  to  50  /xJ  were  focused 
into  the  experimental  cell.  The  pulse  length  was  4  ns, 
and  the  fourth  AS  bandwidth  was  —  1  cm'1. 

The  second  laser  source  of  tunable  193-nm  radiation 
was  a  commercial  two-cavity  excimer  laser  with  a 
grating-tuned  oscillator-amplifier  system  (Lambda  Physik 
EMG150).  The  beam  from  the  oscillator  was  amplified 
with  a  single  pass  through  the  amplifier  cavity  to  an  ener¬ 
gy  of  1-2  mJ/pulse,  with  a  pulse  width  of  12  ns  and  a 
bandwidth  of  0.3  cm-1.  Much  higher  pulse  energies 
could  be  obtained  by  using  unstable  resonator  optics  on 
the  amplifier  and  injection  locking  it  with  the  oscillator 
output,  and  this  configuration  was  used  for  the  two-laser 


LASER  VACUUM  WAVELENGTH  (nm) 


FIG.  3.  Fluorescence  signal  as  a  function  of  ArF  excimer 
wavelength.  These  scans  were  used  to  calibrate  the  wavelength 
scale. 
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depletion  studies  described  below.  The  single-pass 
configuration  eliminated  effects  from  unlocked  broad¬ 
band  output  on  the  two-photon  absorption  measure¬ 
ments,  which  are  sensitive  to  the  laser  bandwidth. 

The  laser  energy  at  193  nm  was  varied  by  admitting  5 
mol  °7c  N:0  in  Ar  into  a  50-cm-long  cell  equipped  with 
Brewster  windows  which  the  laser  beam  passed  through 
prior  to  entering  the  fluorescence  cell.  (5  mol%  Cl2  in 
Ar  was  used  to  attenuate  the  355-nm  beam.)  A  few  hun¬ 
dred  Torr  of  the  N,0  mixture  was  sufficient  to  extinguish 
the  laser  beam,  and  the  attenuation  could  be  continuously 
varied  by  pumping  out  the  mixture.  The  laser  energy  was 
measured  at  the  cell  exit,  with  a  correction  made  for  the 
window  loss  in  the  single  laser  experiments. 

The  wavelength  of  the  excimer  laser  was  calibrated  by 
observing  the  fluorescence  excitation  spectrum  with  both 
the  E,F  (c'  =  6)— X  (u"  =  0)  0(1)  and  Q( 0)  lines  shown  in 
Fig.  1.  Typical  data  for  this  calibration  scan  are  given  in 
Fig.  3.  Additional  wavelength  calibration  experiments 
were  performed  using  1  + 1  two-photon  ionization  in 
NO.14 

Radiation  at  355  nm  was  obtained  from  frequency  tri¬ 
pling  for  the  Nd:YAG  laser.  This  produced  30  mJ  in  a 
pulse  5  ns  long.  The  beam  had  a  roughly  top-hat  spatial 
profile,  but  exhibited  considerable  spatial  mode  structure. 

The  experimental  cell  was  equipped  with  parallel  elec¬ 
trodes  running  the  length  of  the  cell  for  ion  detection, 
and  a  fast  aspheric  lens  to  collect  the  near-ir  fluorescence 
from  the  excited  molecules  due  to  the  transition  E,F 
( L’ '  =  6,  J'=  1 )— »B  (v"  =  0,  J"  =  0,2).  H2  flowed  through 
the  cell  at  pressures  of  0. 2-0. 3  Torr.  We  detected  a  mul¬ 
tiphoton  ionization  signal.  However,  small  quantities  of 
easily  ionized  impurities  in  the  gas  system  (presumably 
vacuum  pump  oil)  interfered  with  the  resonant  ionization 
signal  giving  a  relatively  large  background  signal. 

Fluorescence  collimated  by  the  aspheric  lens  and 
transmitted  by  a  750-nm  bandpass  filter  was  refocused  by 
a  second  lens  (nominal  /  number  —0.7)  onto  an  RCA 
C31034A  photomultiplier  tube.  The  fluorescence  signal 
was  integrated  with  an  Ortec  113  preamplifier  for  all  the 
two-photon  absorption  measurements.  The  fluorescence 
signal  was  recorded  while  scanning  the  laser  through  the 
£,F(u’  =  6)—X  (u"=  0)  (2(1)  transition  at  various  pulse 
energies. 

The  excited-state  decay  rate  was  measured  by  fitting  a 
single  exponential  to  the  time-resolved  fluorescence  signal 
(acquired  on  a  Tektronix  2430  digital  oscilloscope  in 
averaging  mode  without  the  Ortec  113  preamplifier). 
This  was  used  for  the  total  decay  rate  2u"  +  Q  [see 

Eq.  (1)].  The  collisional  quenching  rate  by  H2  was  mea¬ 
sured  to  oe  1.95X  10"’  cm’s'1,  close  to  the  value  previ¬ 
ously  measured  by  Kligler  and  Rhodes. 15 

The  spatial  beam  piofile  was  measured  with  a  pyroelec¬ 
tric  detector  array  which  has  100-/jm  resolution  (Spiricon 
LP-256-1 1-SP).  Scanning  a  pinhole  across  the  beam  gave 
the  same  result  as  the  more  convenient  pyroelectric  ar¬ 
ray.  The  temporal  pulse  profiles  were  measured  by 
reflecting  a  fraction  of  the  beam  onto  a  fast  photodiode 
(100-ps  rise  time,  Instrument  Technology  Limited)  and 
observing  the  output  on  a  1-GHz  oscilloscope  (Tektronix 
7104).  Several  pulses  were  photographed  for  each  run, 
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then  digitized  by  hand.  They  were  subsequently  integrat¬ 
ed  (or  squared  and  integrated)  and  averaged. 

Excited-state  absorption  cross  sections  were  measured 
in  a  pump-probe  arrangement  in  which  the  fourth  AS 
beam  pumped  some  molecules  into  the  E,F  state  and  ei¬ 
ther  the  excimer  beam  or  the  tripled  Nd:YAG  beam  dep¬ 
leted  a  fraction  of  the  excited  population  after  a  several 
nanosecond  time  delay. 13  The  «wo-laser  fluorescence  de¬ 
pletion  experiments  had  the  advantage  that  the  absolute 
number  density  of  excited  molecules  did  not  need  to  be 
known,  since  the  fractional  depletion  of  the  lower  level  in 
the  absorption  of  the  third  photon  was  directly  measured. 
A  short  focal  length  (/=  20  cm)  lens  was  used  in  the 
pump  laser  beam,  generating  a  high  density  of  excited 
molecules  in  a  small  volume  which  was  viewed  by  the 
photomultiplier.  A  significant  fraction  of  the  excited 
molecules  absorbed  a  third  pump  laser  photon.  Because 
the  absorption  cross-section  measurement  did  not  depend 
on  the  excited-state  density,  this  did  not  affect  the  result. 
A  second  laser  pulse  spatially  overlapped  the  pump  laser 
focus,  and  was  delayed  in  time  by  —  10  ns.  [The  excimer 
laser  was  tuned  to  a  wavelength  near  the  E,  F  ( v '  =  6  )•«— X 
(u"  =  0)  Q(  1)  two-photon  transition.  No  fluorescence 
could  be  detected  with  the  pump  beam  blocked.  The 
wavelength  of  the  excimer  laser  beam  was  then  within 
0.05  nm  of  193.26  nm.]  In  order  to  be  certain  that  all  the 
excited  population  was  subjected  to  the  same  ionizing 
laser  fluence,  the  time-delayed  probe  beam  was  not  fo¬ 
cused.  High  pulse  energies  were  then  required  to  pro¬ 
duce  measurable  depletion.  By  placing  a  gate  on  the 
fluorescence  signal  —20  ns  after  the  delayed  laser  pulse 
had  ended,  the  fluorescence  was  observed  to  be  depleted 
by  the  delayed  laser.  A  second  gate  was  placed  at  the 
peak  of  the  fluorescence  signal  prior  to  the  delayed  probe 
laser  pulse;  this  provided  a  signal  proportional  to  the  den¬ 
sity  of  two-photon  excited  molecules  for  shot-to-shot  nor¬ 
malization.  This  signal  was  unaffected  by  the  probe  laser. 
The  ratio  (delayed  fluorescence  signal)/(peak  fluorescence 
signal)  at  each  delayed  probe  laser  power  was  divided  by 
the  value  of  the  ratio  with  the  delayed  probe  laser  absent. 


FIG.  4.  Fluorescence  depletion  measurement  at  193  nm.  The 
solid  curve  shows  a  least -squares  fit  to  an  exponential  decay. 
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FIG.  5.  Fluorescence  depletion  measurement  at  355  nm.  The 
solid  curve  shows  a  least-squares  fit  to  an  exponential  decay. 


Because  the  decay  rate  of  the  excited  population  by 
fluorescence  or  collisional  quenching  was  unaffected  by 
the  probe  laser,  this  normalization  eliminated  effects  of 
those  processes  on  the  results. 

The  pulse  energy  of  the  delayed  laser  was  varied  from 
zero  to  full  power,  and  a  plot  of  delayed  fluorescence 
versus  delayed  laser  was  recorded.  Typical  depletion 
plots  are  shown  in  Fig.  4  for  193-nm  absorption,  and  in 
Fig.  5  for  355-nm  absorption.  The  plots  show  an  ex¬ 
ponential  decay.  The  decay  constant  of  the  delayed 
fluorescence  as  a  function  of  the  fluence  of  the  delayed 
laser  (in  photons/cnr )  equals  the  cross  section  in  cm:  for 
absorption  of  a  photon  from  the  delayed  laser  beam.  The 
fluence  of  the  delayed  laser  was  determined  by  removing 
the  fluorescence  cell  and  placing  a  100-^tm-diam  pinhole 
at  the  waist  of  the  focused  pump  beam,  and  then  measur¬ 
ing  the  energy  of  the  delayed  laser  beam  which  went 
through  the  pinhole.  Since  the  delayed  laser  beam  was 
not  focused,  the  100-/jm  pinhole  had  adequate  resolution 
to  measure  the  laser  fluence. 

B.  Measurements 

Plots  of  Inipeak  signal*  versus  lntlaser  energy)  for  one- 
laser  data  were  prepared  to  determine  the  effective  order 
of  nonlinearity,  k,  where  k  is  the  apparent  number  of 
photons  involved  in  a  given  process: 

k  =3  ln(  signal  )/31n(  energy)  . 

The  plot  of  the  fluorescence  data  had  a  slope  of  2,  which 
is  expected  for  a  two-photon  process  in  the  absence  of 

saturation  Fig.  6  . 

Art  effective  two-photon  excitation  cross  section 
jo  •  ( >  ]  w.:»  calculated  using  the  formula1" 

s  -i-  r  lh  >' 

<7  - - - - -  ,  (1) 

E  -4„,  A  „DJ  F-\t)di 

where  <r  ;  i<j>  (also  known  as  the  two-photon  rate 
coefficient1'’!  is  the  effective  two-photon  cross  section  in 
units  of  cmJs:  SUol/E2  is  the  fluorescence  signal  in  volts 
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FIG.  6.  Intensity  dependence  of  the  fluorescence  signal  for 
data  given  in  run  2  of  Table  II. 


divided  by  the  square  of  the  pulse  energy:  r0  is  the  1  /e~ 
radius  of  the  laser  beam  intensity  (fit  to  a  Gaussian 
profile;  averaged  over  the  fluorescence  viewing  region): 
fui)  is  the  energy  of  a  single  laser  photon,  JV0  is  the  density 
of  molecules  in  the  initial  state  of  the  two-photon  transi¬ 
tion;  and  D  is  the  calibration  constant  for  the  fluores¬ 
cence  detection  system."  Fit )  is  the  measured  temporal 
laser  profile,  defined  as  in  Ref.  12: 

I(r,z,t)=  U(r,z)Flt )  , 

r  Fit)dt=l  , 

w'here  / (r,z,/)  is  the  beam  intensity  and  Uir,z)  is  the 
fluence.  The  fluorescence  quantum  yield  is 

^60 

2  A*"  +  Q  ' 

v" 

where  Aw  is  the  Einstein  A  coefficient  for  the  detected 
fluorescence  in  the  E,F  (v‘  =  6)—*B  (u"  =  0)  band,  Q  is 
the  rate  of  collisional  quenching,  and  £,  •  Abl,~  is  the  sum 
of  Einstein  A  coefficients  for  fluorescence  in  all  bands 
( 2c ■  ^4 6, =  1.01  X  107  sec-1 ).  The  quantity  A^  was  ob¬ 
tained  from  a  theoretical  calculation17  of  the  band  transi¬ 
tion  moment  (Table  I). 

The  fluorescence  detection  system  was  calibrated  with 
spontaneous  Raman  scattering  of  a  visible  laser  beam  by 
a  static  fill  of  ~  1  atm  H2  in  the  experimental  cell  to  pro¬ 
duce  a  known  photon  flux  at  the  wavelength  of  interest. 11 
The  most  important  fluorescence  wavelengths  and  calcu¬ 
lated  branching  ratios  are  listed  in  Table  I.  The  ratio  of 
the  strengths  of  two  rotational  lines  originating  at  the 
same  level  within  a  band  equals  the  ratio  of  the  Hdnl- 
London  factors  for  the  two  lines.  We  collected  photons 
from  two  lines  only,  the  P( 2)  and  R( 0)  lines  originating  at 
the  E,F  (i’’  =  6,  J'=  1)  level  and  terminating  in  the  B 
(r"  =  0,  J"  =  2,0)  levels,  at  wavelengths  752.7  and  746.2 
nm,  respectively.  The  ratio  of  the  calculated  Honl- 
London  factors  for  these  two  lines  in  a  2-2  transition  is 
2. 18  The  calibration  constant  D  [see  Ref.  11,  Eq.  (1)]  was 
measured  in  each  run  at  an  intermediate  wavelength, 
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TABLE  I.  Fluorescence  wavelengths,  Einstein  A  coefficients,  and  branching  ratios  for  H2  E,F 
(u'  =  6, •/’=!). 


Lower  level 

A[R(0>] 

(nm) 

L[/>(2)] 

(nm) 

(s'1)* 

Branching  ratio 

B  <r"=0> 

746.2 

752.7 

3.6X106 

0.36 

B  («.•’*=  1) 

827.6 

835.2 

4.9  X  lO* 

0.49 

B  (i'"  =  2) 

925.7 

934.9 

2.6  X105 

0.03 

B  lc"  =  3> 

1047 

1058 

2.4  X104 

0.002 

B  (u"  =  4) 

1070 

1082 

7.7  XIO4 

0.008 

B  ( v"  =  5> 

1396 

1415 

4.6  X  10- 

0.05 

All  others 

>  1400 

>  1400 

7.3  X  10? 

0.07 

Total 

1.01X10’ 

1.00 

“Reference  17. 


748.8  nm;  then  D  was  measured  once  at  several  wave¬ 
lengths.  Combining  the  ratio  of  line  strengths  and  the 
wavelength  dependence  of  the  apparatus  responsivity,  it 
was  found  that  Dcfr  =  0.S9Dus  g,  where  2>eff  =  f.D752  7 

+  7^746. 2  - 

The  two-photon  excitation  cross  section  was  measured 
several  times,  varying  the  beam  diameter  at  the  focus. 
The  measured  parameters  for  Eq.  (1)  are  listed  in  Table  II 
for  each  run.  In  run  6,  a  cylindrical  beam  profile  was  ob¬ 
tained  by  sending  the  collimated  laser  beam  through  an 
aperture  just  before  the  fluorescence  cell,  in  which  case 
the  radius  of  the  cylindrical  profile  could  be  directly  sub¬ 
stituted  for  r0  in  Eq.  (1).  [See  the  derivation  of  Eq.  11)  in 
Ref.  12.]  Large  variations  between  runs  in  the  photomul¬ 
tiplier  calibration  constant  D  were  caused  primarily  by 
the  use  of  different  photomultiplier  tube  bias  voltages  to 
obtain  optimum  signals  for  data  acquisition. 

The  two-laser  absorption  cross-section  measurement  at 
193  nm  was  performed  several  times  in  a  two-day  period. 
Successive  runs  agreed  to  within  about  2%,  while  all  the 
data  agreed  to  within  10%.  We  attribute  the  10%  spread 
to  random  error  in  calibrating  the  fluence  of  the  delayed 
laser.  An  uncertainty  of  20%  for  the  fluorescence  de¬ 
pletion  measurements  includes  possible  systematic  errors 
(mainly  pyroelectric  detector  calibration).  We  found 


an  excited-state  absorption  cross  section  of 
(6.4+1. 3)X  10“ 18  cm2  at  193  nm. 

The  two-laser  absorption  cross-section  measurement  at 
355  nm  showed  a  greater  run-to-run  variation  (  —  15%) 
than  the  193-nm  measurement.  This  was  probably  due  to 
the  pronounced  spatial  mode  structure  of  the  Nd:YAG 
laser  output  which  makes  it  more  difficult  to  uniformly 
irradiate  the  excited  molecule  region  with  the  355-nm 
laser  output.  We  obtained  a  value  of  (9.7±2.4)X  10“ 18 
cm2  for  the  excited-state  absorption  cross  section  at  355 
nm. 


C.  Role  of  photon  statistics 


The  line-shape-independent,  two-photon  excitation 
cross  section  cr],2  ',  which  is  a  fundamental  property  of  the 
molecule,  is  related  to  the  effective  cross  section  defined 
by  Eq.  (1)  as  follows:16 

dn\<j>)=o$]g(o>)Gn' ,  (2) 


where  g(w)  is  the  area-normalized  line-shape  function, 
and  G!21  is  the  second-order  intensity  autocorrelation 
function  of  the  laser, 


C<2>- 


</2(l)> 

<fU))2  ' 


(3) 


TABLE  II.  Experimental  parameter  values  for  Pq'  measurement. 


Parameter  (units) 

Run  1 

Run  2 

Run  3 

Run  4 

Run  5 

Run  6 

S/E:  i  10'  V/J-’l 

68.4 

7.04 

0.395 

0.217 

0.372 

2.42 

cy,  <10  cm- ) 

1.81 

0.848 

2.03 

1.57 

1.73 

3.39 

-4 «./  |2-4.,  +<2] 

0.114 

0.109 

0.115 

0.114 

0.135 

0.110 

.V„  <  I01'  cm  ')4 

7.15 

7.57 

7.04 

7.08 

5.41 

7.41 

D  <  10  V  cm  sr) 

168 

15.2 

4.94 

1.25 

2.71 

80.5 

/  F:indi  '  10'  s  ') 

3.01 

3.00 

0.755 

0.800 

0.772 

0.871 

<7  <0>  <  10  4  cm4 s) 

1.25 

0.660 

1.11 

1.74 

1.76 

0.599 

G  ’  a  <  10  cm4 1 

4.3 

2.0 

3.6 

3.9 

5.8 

1.3 

Laser  type 

dye 

dye 

ArF 

ArF 

ArF 

ArF 

G  'assumed! 

1.4 

1.4 

2.0 

2.0 

2.0 

2.0 

o  i  10  "  cm4 1 

3.1 

1.4 

1.8 

1.9 

2.9 

0.65 

“Density  in  the  J"~  I  rotational  level  of  H;.  The  quenching  rate  Q  was  calculated  using  the  total  densi¬ 
ty  with  .V,  =  I.52.V,,. 
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where  /(/)  is  the  true  area-normalized  temporal  profile  of 
the  laser  pulse  measured  with  an  infinitely  fast  detector. 
The  angular  brackets  denote  time  averaging  over  a  scale 
long  compared  to  the  coherence  time  of  the  laser,  but 
short  compared  to  the  laser  pulse  length:19 

at  J  i 

</:(r)>  =  -^-  f  ,  +  a'/2(t')dr'  , 
at  J  t 


The  response  time  of  the  photodiode  used  is  in  the  proper 
range;  thus  <  fit))  equals  the  observed  Fit).  The  band¬ 
width  of  the  excimer  laser  was  considerably  smaller  than 
that  of  the  Raman-shifted  laser,  but  the  observed 
linewidths  [and  hence  g(&>)]  were  mainly  determined  by 
the  sizable  Doppler  width  for  H;  (0.90  cm' 1 ).  The  quan¬ 
tity  G  :af  is  obtained  in  absolute  units  by  integrating 
a  :  W)  over  the  excitation  profile,  as  discussed  previous¬ 
ly.1-19 

The  factor  G2'  characterizes  the  photon  statistics  of 
the  pump  beam,  which  must  be  taken  into  account  when 
a  multiphoton  cross  section  is  compared  with  theory. 
Our  detection  apparatus  is  too  slow  to  capture  the  inten¬ 
sity  variations  (on  a  time  scale  of  ps)  of  the  multimode 
laser  output.  The  factor  G‘2  may  be  thought  of  as  ex¬ 
pressing  the  amount  of  “spikiness"  in  the  beam  intensity 
that  cannot  be  resolved  by  the  photodiode.  This  point 
has  been  more  fu  v  discussed  for  two-photon  excitation 
in  atomic  oxygen  in  Ref.  19. 

We  used  G  ;  =  1.4  for  the  Raman-shifted  dyer  laser, 19 
and  2.0  for  the  excimer  laser  (i.e.,  the  excimer  laser  was 
assumed  to  be  a  chaotic  source).  These  choices  are 
reasonably  consistent  with  the  present  results:  the  aver¬ 
age  value  of  G':'<7o21  for  the  Raman-shifted  laser  mea¬ 
surements  (runs  1  and  2)  is  3.1  X  10” 36  cm4,  while  for  the 
excimer  laser  (runs  3-6)  it  is  3.7X10-36  cm4.  The  ratio 
3.7/3. 1  =  1.2,  which  is  fairly  close  to  2.0/1. 4  =1.4.  The 
scatter  is  too  large  to  consider  the  difference  significant. 
The  average  of  G'2'o'q'  from  runs  1  and  2  divided  by  1.4, 
and  runs  3-6  divided  by  2.0,  is  <7o2l  =  (2.0±0.9)X  10-36 


III.  THEORY 


A.  Comparison  of  aj,2’  with  previous  calculations 


The  theoretical  value  of  the  integrated  two-photon  ex¬ 
citation  cross  section  Op2’  (in  cm4)  was  obtained  from  the 
calculated  two-photon  transition  moment3  using  the  for¬ 
mula16 


Oo2l  =  (27T):l 


2 

(W/>2  , 


(4) 


where 

p  _  y  (f  \r\k  )(k\r\g  ) 
*  Eg—  +fui> 


(5) 


is  the  two-photon  transition  moment  for  excitation  by 
one  laser  beam;  here  g,  k ,  and  /  refer  to  the  ground,  vir¬ 
tual  intermediate,  and  final  states  of  the  transition,  re¬ 
spectively;  r  is  the  dot  product  of  the  dipole  moment 
operator  with  the  laser  beam  polarization;  and  at  is  the 
angular  frequency  of  the  laser.  It  is  necessary  to  be  care¬ 
ful  in  applying  the  results  of  Ref.  5.  When  the  value  of 
P 2  from  Ref.  5  of  19.7  a.u.  is  entered  into  Eq.  (4),  the  re¬ 
sult  is  much  larger  than  our  measured  value.  However, 
the  authors  of  Ref.  5  originally  calculated  a  transition 
moment  Af  for  absorption  of  one  photon  from  each  of 
two  statistically  independent  laser  beams: 


(f\r\k  )  ( k\r'\g  ) 
Eg  —  Ek  —  fiat 


<  f\r’\k  )  <  /cl rjg  ) 
Eg  —  Ek  —  flat 


(6) 


(where  the  prime  denotes  the  second  laser  beam).  They 
then  set  r'  =  r  and  at'=at.  To  avoid  double  counting,  the 
resulting  transition  moment  must  be  divided  by  2  in  or¬ 
der  to  use  it  in  Eq.  (4). 20  Since  the  transition  moment  is 
squared  in  Eq.  (4),  the  cross  section  is  reduced  by  a  factor 
of  4.  The  correctly  calculated  cross  section,  Op2 
=  2.8X10~36  cm4,  is  in  reasonable  agreement  with  our 
measurements. 


B.  Direct  photoionization  of  the  E, Estate 

Several  processes  fall  under  the  general  category  of 
photoabsorption  by  the  E,F  state.  In  this  section  we  con- 


TABLE  Ill.  Calculated  cross  sections  for  H:  (E.F,  v=6)+fja>. 


Process 

<7j5.  (10  "  cm2) 

<7|»,  ( 10  18  cm2) 

Direct  photoionization 

to  l’  *.  =  1 

0.37 

0.10 

to  u  .  =  2 

8.11 

2.03 

to  t'-  =13 

0.52 

Sum  over  all 

open  channels 

8.50 

3.15 

Photoexcitation  of  '2„ 

autoiomzation 

0 

7.13 

dissociation  (neutral-)- ionic) 

3.55 

0.52 
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sider  direct  (i.e.,  conventional)  photoionization: 

H:*+fcu— H2++e~  , 

which  is  described  by  the  matrix  element  between  the 
bound  E,F  state  and  the  electronic  continuum  plus  the 
bound  H2  +  state.  The  cross  section  for  this  process  may 
be  estimated  from  calculations  in  the  literature.  The 
work  of  Cohn21  provides  tabulated  cross  sections  as  a 
function  of  the  ejected  electron  energy.  We  interpolated 
these  and  multiplied  by  the  appropriate  Franck-Condon 
factors.  The  results  are  tabulated  in  Table  III.  This  ap¬ 
proximation  neglects  the  R  dependence  and  the  electron 
energy  dependence  of  the  transition  moment. 

More  sophisticated  calculations  of  the  direct  photoion¬ 
ization  process  have  been  carried  out  by  Rudolph  et  al. 22 
These  authors  reported  relative  cross  sections  for  the  final 
vibrational  states,  but  not  the  absolute  numbers.  An  esti¬ 
mate  of  the  direct  photoionization  cross  section  from  the 
E,F  (u  =  6)  state  to  H,*(u=2)  obtained  from  their 
work20  is  1.5X10-18  cm2,  which  is  about  25%  smaller 
than  the  corresponding  value  in  Table  III. 

C.  Photoexcitation  of  the  ’2.  autoionizing  state 
from  the  E,  F  state 

Since  the  measured  cross  section  for  absorption  by  the 
E,F  lr'  =  6)  state  is  larger  than  the  calculated  upper  limit 
for  direct  molecular  photoionization,  another  mechanism 
may  play  a  role.  Recent  work  in  this  laboratory9,23  and 
elsewhere24  2?  has  demonstrated  that  excitation  to  a  dou¬ 


bly  excited,  dissociating  autoionizing  state  is  an  impor¬ 
tant  channel  in  the  photoionization  of  the  C  1riB  state  of 
H2.  We  therefore  performed  quantitative  calculations  to 
assess  the  role  of  this  channel  in  the  present  situation. 

The  mechanism  considered  is  illustrated  in  Fig.  7.  The 
initial  level  is  u'  =  6  of  the  E,F  state.  Excitation  to  the 
lowest  ’2*  autoionizing  stele  is  shown  in  the  figure.  This 
state  has  been  calculated  by  Guberman26  and  its  electron¬ 
ic  configuration  is  2pau2sag.  It  is  the  most  likely  candi¬ 
date  final  state,  because  it  is  reached  by  a  strong  one- 
electron  transition  2pcru*-\sog  from  the  inner-well  E 
lsog2sog  component  of  the  E,F  state,  and  by  a  one- 
electron  Iscjg*—  2pau  transition  from  the  outer-well  F 
(2 pa  u  >2  component  of  the  E,F  state.  Once  excitation  has 
occurred,  the  nuclei  start  to  move  outward  on  the  '2* 
potential  curve.  This  electronic  state  is  unstable  with 
respect  to  autoionization,  which  occurs  as  the  nuclei  dis¬ 
sociate.  This  process,  electron  emission,  leads  to  the  final 
state  H2+  +e  If  autoionization  does  not  occur,  the  nu¬ 
clei  continue  to  follow  the  '2^  curve  and  ultimately 
reach  an  asymptotic  state  of  the  form  H(ls)  +  H*(n/),  or 
H+  +H~. 

The  process  just  described,  which  leads  to  molecular 
autoionization  IH/+f),  neutral  dissociation  (H*+H). 
or  ionic  dissociation  (H  *  +  H "  ),  was  treated  using  the 
technique  developed  by  Hickman.9  The  method  treats 
the  dynamics  of  nuclear  motion  by  solving  Schrodinger’s 
equation  with  a  nonlocal,  complex  potential.  Cross  sec¬ 
tions  for  total  dissociation  and  molecular  autoionization 
were  obtained,  as  well  as  the  distribution  of  final  vibra- 


CROSS  SECTION  (to-’8  cm2)  R  (units  of  a<j) 

FIG.  7.  (a)  Calculations  for  direct  photoionization  ( . ),  photodissociation  through  the  '2„  state  ( - ),  and  autoionization 

i - l  through  the  '2„  state  of  H;.  (b)  Potential  curves  of  Hj  illustrating  the  mechanisms  contributing  to  the  calculations  in  the  left 

panel.  Notice  that  the  energy  (vertical)  scales  are  the  same  in  both  panels,  thereby  clearly  exhibiting  the  threshold  behavior  of  all  the 
cross  sections. 
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tional  states  resulting  from  molecular  autoionization. 
The  necessary  quantities  for  the  calculation  are  the  fol¬ 
lowing:  (1)  the  transition  matrix  element  p(R)  connect¬ 
ing  the  initial  E,F  state  with  the  '2*  state  ( R  is  the  inter- 
nuclear  distance);  (2)  the  potential  curve  for  the  '2* 
state;  and  (3)  the  matrix  element  Fd(R)  connecting  the 
'2J  state  with  the  final-state  continuum  H2++e_.  The 
determination  of  these  potentials  and  matrix  elements  is 
described  in  the  following  paragraph. 

To  calculate  the  transition  moment  /x(/?)  connecting 
the  '2*  autoionizing  state,  one  must  take  into  account 
the  R  dependence  of  the  E,F  electronic  wave  function. 
This  can  be  done  in  an  approximate  way  by  considering 
the  two  dominant  electronic  configurations,  lsag2sog 
and  F2pa u  )2.  Each  of  these  configurations  corresponds  to 
a  diabatic  potential  curve;  the  avoided  crossing  of  these 
curves  leads  to  the  well-known  double-well  nature  of  the 
E,F  state.  We  can  approximate  the  adiabatic  electronic 
wave  function  as  follows: 

f(R  )  =  (cos0t/?)](  lsog2sog  )  +  [sin0(/l!  )](2po„  )2  . 

(7) 

Estimates  of  the  mixing  angle  0  may  be  obtained  from 
the  diabatic  potential  curves  and  coupling  matrix  ele¬ 
ments  determined  by  Hazi  et  al.21  and  Ross  and 
Jungen.2*  One  sets  up  a  two-state  Hamiltonian  in  the 
standard  way:  Hlt(R )  is  the  energy  of  the  (2 pau  )2  state; 
H::  is  the  energy  of  the  lsog2sag  state;  and 

//i:  =  i-0/[2-6s(R)]v:  , 

where  r(l  =  0.048  a.u.,r  and  &,(/?)  is  the  quantum  defect 
function.  There  is  some  uncertainty  in  the  estimate  of 
H  ]2  because  it  was  necessary  to  extrapolate  u0  beyond  the 
range  in  which  it  was  calculated  by  Hazi  et  al.  The  re¬ 
sult  is 

tan[20(Ri]  =  2//i:/(//u-tf22)  . 

We  must  now  evaluate 

p[R)=('\f E  f  {Ryp:.2po u2so g ) 

=  [cos0(R )](  \sa g\r  cosd\2pa u ) 

+  [sin©(R)](2pali!rcos0i2sa8)  ,  (8) 

where  r  and  6  are  spherical  electronic  coordinates. 

The  first  matrix  element  is  for  the  2po„*—  ls(7g  core 
transition  and  has  been  evaluated  by  Bates. ‘9  Its  value  is 
—  R/ 2  times  a  correction  factor  close  to  unity.  The 
second  matrix  element  is  for  the  2sog—  2pou  transition. 
It  was  evaluated  numerically  using  single  center  wave 
functions  determined  from  quantum  defects.  The  results 
are  well  approximated  by  the  polynomial  1.571 
-0. 3807 R  +0.062  68 R  ;  (atomic  units  are  usedi. 

The  '2*  potential  curve  has  been  calculated  by  Guber- 
man.a  The  matrix  element  Fd(R)  is  obtained  from  the 
resonance  width  HR)  calculated  by  Tennyson  and  No- 
Mv,  '  using  the  relation 

rrl(R)=-[r(/?)/2r]‘  2  . 


Tennyson  and  Noble  calculated  HR )  for  R  <  2.8  a.u.;  for 
larger  values  we  took  the  value  at  R  =  2.8  a.u. 

With  these  potential  curves,  and  the  estimates  of  the 
matrix  elements,  we  can  treat  photoexcitation  of  the  1 2„ 
autoionizing  state.  The  calculated  results  are  presented 
and  compared  with  our  experimental  measurements  in 
Sec.  IV. 

IV.  DISCUSSION 

The  measurements  of  the  two-photon  excitation 
( E,F—X )  cross  section  are  in  reasonable  agreement  with 
theory,  as  discussed  in  Sec.  Ill  A.  The  excitation  process 
is  straightforward  and  appears  to  need  no  further  discus¬ 
sion. 

Photoabsorption  by  the  E,F  (v'  =  6)  state  is  more  com¬ 
plex,  because  several  alternate  channels  are  available.  We 
write  several  possibilities  schematically  as  follows: 


H^+e-  ,  (9a) 

H**— ►H2t  +  e ~  ,  (9b) 

H2*(£,F)  +  *G>—  H:*— H*  +  H  ,  (9b’) 

HJ*— ‘H*  +H”  ,  (9b”) 

H+H++e~  .  (9c) 


Some  of  these  channels  have  already  been  observed  or 
discussed  in  the  literature.  Direct  photoionization,  (9a), 
is  well  known.22  Channels  (9b),  (9br),  and  (9b")  denote 
the  excitation  and  subsequent  decay  of  highly  excited  au¬ 
toionizing  states.  Competition  between  autoionization 
and  dissociation  of  these  channels  has  been  studied.  '1 
Channel  (9c)  is  possible  for  sufficiently  energetic  photons, 
and  has  been  inferred  from  observations  of  proton  pro¬ 
duction.  10  A  complicating  feature  of  the  experimental 
analysis  of  this  channel  is  that  other  possible  sources  of 
H  +  ions  are  photoionization  of  the  excited  atomic  hydro¬ 
gen,  H*,  produced  by  (9b'),  or  photodissociation  of  H2*. 

Previous  experiments  have  provided  relative  cross  sec¬ 
tions  for  several  of  the  channels  listed  in  (9).  For  exam¬ 
ple,  the  photoelectron  energy  distribution  has  been  mea¬ 
sured  for  (9)  at  193  nm  and  provides  evidence  of  a  strong¬ 
ly  non-Franck-Condon  distribution  of  final  vibrational 
states  of  H2+. 32,33  These  results  suggest  that  (9a)  is  not 
the  only  important  channel.  The  present  work  provides 
complementary  information,  because  absolute  cross  sec¬ 
tions  are  measured  and  compared  with  theoretical  calcu¬ 
lations.  Since  the  measured  cross  section  for  photoab¬ 
sorption  was  significantly  larger  than  the  theoretically 
calculated  cross  section  for  channel  (9a),  additional  chan¬ 
nels  probably  play  a  role. 

Theoretical  calculations  were  carried  out  for  compar¬ 
ison  with  the  experimental  data,  which  were  obtained  at 
values  of  fui)  corresponding  to  355  and  193  nm.  The  cal¬ 
culations  treat  channels  (9a),  (9b),  (9b’l,  and  (9b”l,  for 
which  estimates  of  the  coupling  terms  could  be  obtained 
as  described  in  Sec.  III.  These  calculations  suffice  to 
demonstrate  the  importance  of  channels  (9b),  (9b  ),  and 
(9b”).  Coupling  terms  to  treat  channel  (9c)  have  not  been 
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TABLE  IV.  Comparison  of  experimental  and  theoretical  cross  sections. 


a\(E,F  (i>'  =  6)«— (u"  =  0)) 

<7jm(£,F  (u"  =  6)-t- 193  nm) 

ottn{E,F  (u''  =  6)  +  355  nm) 

HO" 36  cm*) 

( 10“ 18  cm2) 

(10" 18  cm2) 

Experimental 

2.0±0.9 

6.4±  1.3 

9.7±2.4 

Calculated 

2.8 

10.8 

12 

calculated,  and  thus  only  qualitative  comments  can  be 
made  about  this  channel. 

The  theoretical  calculations  are  summarized  in  Table 
III  and  in  Fig.  7.  An  interesting  implication  of  the  calcu¬ 
lations  is  that  the  relative  importance  of  the  various 
channels  (9)  changes  dramatically  with  photon  energy. 
This  behavior  can  be  best  understood  by  examining  Fig. 
7.  Cross  sections  for  the  various  channels  are  shown. 
The  short-dashed  curve  is  direct  photoionization,  and  is 
summed  over  all  open  final  vibrational  channels.  It  is 
shown  for  energies  above  the  threshold  for  excitation  of 
the  i'  =  2  vibrational  state  of  H2*.  The  jump  at  about 
40000  cm"1  corresponds  to  the  opening  of  the  u=13 
state  of  the  ton.  The  Franck-Condon  factor  for  this  state 
(0.089)  is  second  only  to  that  for  the  v=2  state  (0.754), 
and  comes  from  the  portion  of  the  wave  function  in  the 
outer  F  well.  The  long-dashed  curve  is  photodissociation 
(including  both  neutral  and  ionic  channels)  on  the  '2,, 
autoionizing  state.  The  details  of  this  curve  are  uncertain 
because  of  the  approximate  transition  matrix  element 
used,  but  the  general  features  make  physical  sense.  The 
threshold  occurs  at  the  energy  at  which  H(n  =2)  +  H 
may  be  formed.  The  oscillatory  structure  reflects  the  os¬ 
cillations  of  the  vibrational  wave  function  in  the  outer 
well.  In  the  region  of  oscillations,  photodissociation  is 
effectively  the  only  available  channel.  The  overlap  be¬ 
tween  the  dissociating  wave  function  (H*  +  H)  and  the 
H:"  vibrational  wave  functions  is  negligible  up  to  about 
30000  cm"1.  Above  this  energy,  the  solid  curve  starts  to 
rise.  This  curve  corresponds  to  photoexcitation  of  the 
'2U  state,  followed  by  autoionization,  which  can  only  be 
large  when  the  wave  function  of  the  dissociating  state 
(determined  by  the  photon  energy)  has  a  good  overlap 
with  the  wave  functions  of  the  open  vibrational  channels. 
It  ts  interesting  to  note  that  the  oscillations  in  the  photo- 
dissociation  curve  are  quenched  once  autoionization  be¬ 
comes  probable. 

These  calculations  suggest  the  following  interpretation 
of  the  observed  photoabsorption  signal.  At  355  nm,  the 
process  is  predominantly  direct  photoionization,  with  a 
small  component  of  photodissociation.  Dissociative  ion¬ 
ization.  <9ci.  and  ionic  dissociation,  (9b”),  are  not  ener¬ 
getically  allowed.  At  193  nm,  direct  photoionization  is 
much  smaller,  and  photoexcitation  of  the  state  fol¬ 
lowed  by  autoiomzation  is  quite  important.  Dissociative 
ionization  is  possible  at  this  energy,  but  our  present  data 
do  not  distinguish  between  the  production  of  H;*  and 
H  ions  Anderson  cl  al.':  argued  that  less  than  5 %  of 
all  the  ions  were  H",  although  Buck  er  at.10  measured 
this  fraction  to  be  about  33^ .  Further  studies  are  in  pro¬ 
gress  to  investigate  this  point. 

We  find  that  the  sum  of  the  calculated  cross  sections 
for  (9a),  (9b),  (9b  ),  and  (9b”)  compares  reasonably  well 


with  the  measured  data,  as  summarized  in  Table  IV. 
Channels  (9a)  and  (9b!  may  interfere,  so  that  the  total 
cross  section  may  not  be  equal  to  the  sum  of  the  cross 
sections  calculated  for  individual  channels.  We  have  not 
taken  that  into  account.  The  crucial  point  we  wish  to 
make  is  that  the  measured  cross  section  for  photoabsorp¬ 
tion  is  larger  than  the  cross  section  calculated  for  direct 
photoionization  (a  factor  of  2  larger  at  193  nm).  The  ad¬ 
ditional  channels  considered  clearly  provide  a  substantial 
contribution  to  the  total  signal. 

V.  CONCLUSIONS 

The  measured  integrated  two-photon  excitation  cross 
section,  oj/'  =  (2.0±0.9)X  10~36  cm4,  was  in  good  agree¬ 
ment  with  the  calculated  value,  o02,  =  2.8X  10" 36  cm"4, 
for  the  H,  E,F'lg  ( v '  =  6 )—X  1 2g  (u"=0)  Q(l)  transi¬ 
tion.  The  cross  section  for  absorption  by  the  excited  E,F 
state  of  a  third  193-nm  photon,  oabs~(6A-  I.3)X  10" 18 
cm2,  was  found  to  be  substantially  larger  than  our  es¬ 
timated  value  for  direct  photoionization,  3.2  X  10“ 18  cm2. 
Excitation  of  the  '2„  dissociating,  autoionizing  state  was 
treated  and  shown  to  contribute  significantly  to  the  pho¬ 
toabsorption.  Agreement  was  good  between  theory  and 
experiment  for  absorption  by  the  E,F  state  of  a  355-nm 
photon.  The  sum  of  cross  sections  calculated  for  the 
separate  channels  was  12X10“ 18  cm2;  we  measured 
olbs  =  (9.7±2.4)X  10“ 18  cm2. 

The  technique  of  two-photon  excited  fluorescence 
(TPEF!  appears  promising  as  an  analytical  technique  to 
measure  densities  of  vibrationally  excited  H2.  Scaling 
from  our  demonstrated  signal-to-noise  ratios  for  t  "  =  0, 
we  estimate  that  densities  of  1010  cm  “ 3  per  quantum  state 
could  be  detected  by  using  TPEF,  using  a  laser  pulse  of 
10  mJ  focused  to  10“2  cm“2  with  a  10-ns  pulse  length.  If 
the  two-photon  absorption  cross  sections  for  r”  =  1  —  14 
are  as  large  as  the  cross  section  reported  here  for  c"  =  0, 
it  should  be  possible  to  determine  the  vibrational  state 
distributions  in  the  multicusp  H“  ion  sources  currently 
under  development.6  34 

Note  added.  The  theoretical  value  of  <jJ  has  recently- 
been  reevaluated  by  Huo.36  The  new  value  is  3.6  x  10  16 
cm4.  This  will  be  discussed  further  in  subsequent  publi¬ 
cations. 
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Photoelectron  energy  spectra  following  (2  + 1  )-multiphoton  ionization  of  H2  via  the  double-well 
E.F'i*  state  have  been  measured  using  a  magnetic-bottle  electron  spectrometer.  The  range  of  vi¬ 
brational  levels  studied,  uE  F=3-9,  includes  states  localized  at  the  bottom  of  the  E  well,  and  also 
states  that  span  both  wells.  The  branching  among  the  vibrational  ionization  channels  is  governed 
by  the  strong  R  dependence  of  the  electronic  wave  function  and  the  degree  of  localization  of  the  vi¬ 
brational  wave  function  in  the  E  or  the  F  well.  Theoretical  analysis  confirms  that  at  least  two  mech¬ 
anisms  contribute  to  the  observations:  a  direct  process  involving  only  the  E,F  state  and  the  H*  ion¬ 
ization  continuum,  and  an  indirect  process  involving  the  '2*  autoionizing  state.  Mass  spectroscop¬ 
ic  measurements  carried  out  for  vE  f— 9,  J  —  1  show  that  about  25%  of  the  ions  produced  were  H~ . 
Qualitative  arguments  suggest  that  most  of  these  protons  arise  from  dissociative  ionization  to  the 
H2  continuum. 


I.  INTRODUCTION 

Dynamical  processes  in  diatomic  molecules  can  be  in¬ 
vestigated  by  multiphoton  ionization  and  measurement  of 
the  photoelectron  energy.  Such  spectra  directly  provide 
final-state  distributions,  from  which  many  features  of  the 
molecular  interactions  may  be  inferred.  For  example,  re¬ 
cent  studies  of  the  multiphoton  ionization  of  H2  (Refs. 
1-9)  revealed  anomalous  distributions  of  final  vibrational 
states  of  Hi".  Although  surprising  at  first,  these  distri¬ 
butions  were  shown  to  arise1011  from  an  indirect  process 
involving  dissociating,  autoionizing  states  in  the  H2  con¬ 
tinuum.  A  feature  of  the  indirect  ionization  process  is 
that  final  vibrational  states  may  be  populated  that  have  a 
poor  Franck-Condon  factor  with  the  vibrational  level  of 
the  neutral  intermediate  state. 

This  paper  reports  the  extension  of  our  previous  mea¬ 
surements  of  multiphoton  ionization  (MPI)  of  H2  to  the 
case  where  specific  vibrational  levels  of  the  E,F  state  are 
prepared  by  two-photon  excitation  and  are  photoionized 
by  a  third  photon.  The  E,F  states  present  new  features 
because  of  the  unusual  double-well  structure  of  the  poten¬ 
tial  curve,  and  have  been  studied  by  several 
groups.’-5 12-14  The  electronic  structure  of  these  states 
exhibits  a  strong  dependence  on  the  internuclear  separa¬ 
tion  R.  Depending  on  whether  the  vibrational  state  is  lo¬ 
calized  mostly  in  the  E  well  or  mostly  in  the  F  well,  one 
expects  different  channr’s  to  come  into  play.  By  tuning 
the  laser  to  alternate  vibrational  levels  in  the  E,F  state 
i'£  f  we  have  observed  dramatic  changes  in  the  final-state 
distributions.  Furthermore,  theoretical  analysis  suggests 
a  consistent  interpretation  of  this  behavior  in  terms  of  the 
interplay  of  several  alternate  ionization  channels. 

This  paper  is  organized  as  follows.  Section  II  describes 
the  experimental  procedure.  Section  III  reviews  the  in¬ 
direct  autoionizing  channel  and  describes  theoretical  cal¬ 
culations  based  on  the  nonlocal  complex  potential  model 
that  include  this  channel.  Section  IV  describes  and 
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discusses  our  results  and  Sec.  V  presents  concluding  re¬ 
marks. 

II.  EXPERIMENTAL  PROCEDURE 

The  multiphoton  ionization  scheme  is  shown  in  Fig.  1: 
two  photons  in  the  energy  range  around  6.5  eV  resonant¬ 
ly  excite  intermediate  E,F  levels  of  known  rovibrational 
quantum  number.  The  E,F  double-minimum  potential 
contains  low  vibrational  states  which  are  trapped  mostly 
in  either  the  inner  E  well  or  the  outer  F  well,  as  well  as 
vibrational  states  with  their  wave  functions  spread  over 
the  potential  barrier  and  out  to  large  R.  The  third  pho¬ 
ton  from  the  same  laser  ionizes  the  molecule.  Photoelec¬ 
trons  from  the  various  resulting  ionic  states  are  collected 
in  a  time-of-flight  (TOF)  electron  spectrometer. 

We  use  a  YAG  (yttrium  aluminum  garnet)  pulsed  laser 
system  and  a  high-pressure  H2  Raman  cell  to  generate 
tunable  light  in  the  wavelength  range  from  180  to  200  nm 
at  100  /pulse  energy  level.  The  laser  beam  is  focused 
by  a  15-cm  focal  length  lens  and  forms  a  500-jum  size  spot 
at  the  center  of  the  electron  spectrometer.  Simultaneous¬ 
ly  injected  into  the  same  spot  is  a  supersonic  jet  of  H: 
from  a  pulsed  valve.  The  spectrometer  is  of  magnetic- 
bottle-type  design,  energized  by  permanent  magnets.  It 
has  a  2ir  collection  efficiency  due  to  its  angular  integra¬ 
tion  and  a  high  time  and  energy  resolution  due  to  its  field 
design.  The  principle  and  details  of  this  instrument  have 
been  published  by  Tsuboi  et  al. 15  The  spectrometer  is 
differentially  pumped  with  two  turbo  molecular  pumps  to 
10-8  torr,  and  the  average  background  pressure  of  H2  gas 
in  the  magnet  chamber  is  kept  in  the  10-8-torr  range 
during  the  experiment.  Photoelectrons  created  by  laser- 
H2  interaction  follow  the  magnetic  field  lines  into  the  65- 
cm  TOF  tube  and  are  collected  on  a  multichannel  plate 
detector.  We  record  the  photoelectron  signal  with  an 
eight-channel  real-time  time-to-digital  converter  (TDC) 
(Le  Croy’s  model  4208),  which  works  in  a  multihit  mode 
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INTERNJCLEAR  DISTANCE  (A) 

FIG.  1.  Diagram  of  the  relevant  potential  curves  for  the  mul- 
tiphoton  ionization  of  H,  through  the  E,F  state.  Solid  curves 
are  H:  states,  dashed  curves  are  Hi*. 

in  this  experiment.  Its  intrinsic  time  resolution  is  1  nsec, 
which  should  enable  us,  in  principle,  to  resolve  rotational 
splitting  for  4-eV  electrons.  Data  acquisition  is  con¬ 
trolled  by  a  PDP-11  computer.  The  use  of  the  TDC  in 
replacement  of  the  analog  averaging  method  that  was  ap¬ 
plied  in  our  previous6  photoelectron  experiment  involv¬ 
ing  H;  has  greatly  improved  our  signal-to-noise  ratio  and 
time-energy  resolution.  Nevertheless,  we  still  have  to 
subtract  each  data  spectrum  with  its  corresponding  back¬ 
ground  spectrum  taken  when  the  laser  is  detuned  from 
the  resonance  frequency,  owing  to  an  unidentified  impuri¬ 
ty  that  generates  photoelectrons  in  the  energy  range  from 
0.5  to  2  eV. 

Currently  we  do  not  apply  any  external  field  to  extract 
or  retard  electrons.  However,  an  unidentified  source  of 
contact  potential  limits  collection  of  photoelectrons  to 
those  above  400  meV.  The  retardation  of  electrons 
caused  by  the  contact  potential  is  corrected  by  calibrating 
the  spectra  against  the  known  ionization  potential  of  H2 
and  the  vibrational  spacings  of  H;*.  We  have  also  set  up 
a  time-of-flight  ion  spectrometer  to  monitor  the  creation 
of  ion  species  in  the  same  multiphoton  processes. 

III.  THEORY 

Recent  theoretical1011  16  and  experimental6-9  studies 
have  led  to  a  realization  of  the  importance  of  dissociat- 
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ing,  autoionizing  states  in  the  photoionization  of  elec¬ 
tronically  excited  molecules.  Photoexcitation  of  these 
states,  followed  by  autoionization,  provides  an  additional 
contribution  to  photoionization  and  may  lead  to  strongly 
non-Franck-Condon  distributions  of  final  vibrational 
states.  Photodissociation  to  excited  neutral  products  is 
also  possible.  This  section  summarizes  the  mechanism 
responsible  for  these  effects. 

Figure  1  shows  the  relevant  potential  curves  in  the 
present  situation.  We  use  the  term  direct  photoionization 
to  denote  the  one-electron  transition  between  the  initial 
state  H2  E,F  and  the  electronic  continuum  e+H2+.  This 
transition  is  caused  by  the  dipole-matrix  element  cou¬ 
pling  the  excited-state  orbital  and  the  Coulomb  scattering 
function.  An  additional  process  that  must  be  considered 
is  the  photoexcitation  of  the  dissociating,  autoionizing 
*X*  state,  which  has  the  electronic  configuration 
2sag2pau,  and  is  also  coupled  to  the  E,F  state  by  a  one- 
electron  transition.  If  the  transition  to  the  'X*  state 
occurs,  then  the  nuclei  start  to  dissociate  just  as  in  photo¬ 
dissociation.  However,  the  autoionizing  state  is  unstable 
with  respect  to  electron  emission.  As  the  nuclei  move, 
there  is  a  competition  between  autoionization  and  disso¬ 
ciation.  Roughly  speaking,  if  autoionization  occurs  be¬ 
fore  the  nuclei  separate  to  the  R  *  for  which  the  'X*  state 
is  lower  in  energy  than  the  2?.g  state  of  H2+  (~4a0), 
then  the  final  state  will  be  H2+  plus  a  continuum  elec¬ 
tron,  corresponding  to  indirect  photoionization.  For  the 
present  system  and  for  others  studied  in  our  laborato¬ 
ry, 6-11  our  calculations  indicate  that  this  outcome  occurs 
with  at  least  90-95  %  probability.  If  the  nuclei  reach  R  * 
without  autoionization,  they  encounter  a  series  of  curve 
crossings  with  the  Rydberg  potential  curves  converging 
to  the  state  of  H2+.  Depending  on  the  outcome  of 
the  curve  crossings,  the  nuclei  may  ultimately  dissociate 
or  undergo  a  transition  to  a  bound  state  and  vibrationally 
autoionize  or  radiate. 

For  the  present  situation,  the  lifetime  for  electron 
emission  is  on  the  order  of  8  X  10“ 16  sec.17  A  typical  dis¬ 
sociation  time  estimated  by  classical  mechanics  is  about 
4X10-15  sec,  or  about  five  times  longer.  Therefore 
significant  nuclear  motion  is  likely  during  the  time  it 
takes  to  autoionize.  During  this  time,  electronic  energy 
may  be  converted  to  nuclear  kinetic  energy.  Contrast 
this  situation  to  direct  photoionization  in  which  the  elec¬ 
tron  is  emitted  before  the  nuclear  position  and  momen¬ 
tum  have  time  to  change  substantially.  From  this  point 
of  view,  one  can  see  why  the  Franck-Condon  approxima¬ 
tion  is  applicable  to  direct  photoionization,  but  not  to  the 
indirect  process. 

The  mechanism  just  described  has  been  treated  by  solv¬ 
ing  a  modified  Schrodinger  equation  with  a  nonlocal 
complex  potential.  The  first  application"  of  this  method 
was  to  the  photoionization  of  the  C  'llu  state,  and  a  cal¬ 
culation  for  the  photoionization  of  the  H2  (v£  F  =  6)  has 
been  compared  with  other  experiments  from  this  labora¬ 
tory.18  In  this  section,  we  summarize  the  determination 
of  the  necessary  potential  curves  and  describe  calcula¬ 
tions  of  the  distribution  of  final  vibrational  states  for 
several  initial  vibrational  states  of  the  E,F  potential. 
Cross  sections  for  the  direct  photoionization  process  are 
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estimated  and  added  to  those  for  the  autoionization  pro¬ 
cess.  The  results  presented  in  Sec.  Ill  will  show  each 
contribution  independently.  In  principle,  the  processes  of 
direct  photoionization  and  excitation-autoionization  in¬ 
terfere.  One  should  add  complex  amplitudes  and  then 
square,  rather  than  add  cross  sections  as  we  have  done. 
Dixit  et  at. 16  have  carried  out  such  calculations  on  relat¬ 
ed  systems.  However,  because  we  are  invoking  approxi¬ 
mate  coupling  terms  from  several  sources  with  various 
phase  conventions,  attempting  to  treat  interference  does 
not  seem  warranted. 

The  potentials  and  coupling  terms  used  to  treat  the 
phoioexcitation  from  the  £,£  state  are  discussed  in  detail 
elsewhere.  We  briefly  list  these  terms  and  how  they  were 
obtained,  (i)  The  autoionizing  potential  was  calcu¬ 
lated  by  Guberman.19  (ii)  The  width  HR),  or  electron 
emission  rate,  was  taken  from  calculations  of  Tennyson 
and  Noble.17  (iii)  The  transition  moment  p(R)  connect¬ 
ing  the  £,£  state  with  the  state  was  estimated  by 
considering  the  two  dominant  electronic  configurations 
of  the  E,F  states,  as  described  in  Ref.  18. 

The  direct  photoionization  cross  sections  may  be  deter¬ 
mined  from  previous  work.  Cohn20  tabulated  values  for 
several  excited  states  of  H,  using  the  Coulomb  approxi¬ 
mation  for  ihe  continuum  wave  function.  More  recently, 
Rudolph  et  air'  have  carried  out  more  sophisticated  cal- 


tional  states  following  direct  photoionization,  but  they 
did  not  report  the  absolute  values  of  the  cross  sections. 
Their  results  appear  to  indicate22  that  the  calculations  of 
Cohn  are  accurate  to  within  30%.  For  the  present  calcu¬ 
lations,  we  have  evaluated  direct  photoionization  cross 
sections  by  interpolating  Cohn’s  photoionization  cross 
sections  as  a  function  of  photoelectron  energy  and  multi¬ 
plying  them  by  the  appropriate  Franck-Condon  factor. 

In  summary,  we  calculated  the  distribution  of  final  vi¬ 
brational  states  by  the  following  procedure.  We  calculate 
the  cross  section  for  the  direct  photoionization  process, 
using  the  Franck-Condon  approximation.  To  this  we  add 
the  cross  section  for  reaching  the  same  final  state  by  go¬ 
ing  through  the  '2^  autoionizing  state,  using  the  nonlo¬ 
cal  complex  p; ;  r—ntial  method.11  It  is  this  second  contri¬ 
bution  that  causes  the  strongly  non-Franck-Condon  na¬ 
ture  of  the  distribution. 

IV.  RESULTS  AND  DISCUSSION 
A.  Available  channels 

We  have  measured  the  TOF  spectra  of  photoelectrons 
obtained  in  the  multiphoton  ionization  of  H2  through 
several  vibrational  states  of  the  E,F'2*  double-well  po¬ 
tential.  We  first  indicate  schematically  the  open  channels 


culations.  They  calculated  the  distribution  of  final  vibra-  that  could  lead  to  the  electrons: 

| 

H2+(e  *  )  +  e 

H*  +H 

(la) 

H;  E.F‘2.;  vEF,J  =  \)^  1 

h2**('i;)  -■ 

H2+(e  +  )  +  e 

H  +  H+  -t-e 

(lb) 

1 

H+H~+e  . 

(lc) 

W’e  focus  our  attention  on  the  channels  available  after  the 
first  two  photons  have  been  absorbed  to  form  the  E,F 
state.  Channels  (lal  and  (lb),  direct  photoionization,  and 
excitation  of  the  autoionizing  state,  respectively,  were 
discussed  in  Sec.  III.  Channel  (lc)  is  dissociative  ioniza¬ 
tion,  and  will  be  discussed  in  more  detail  below. 

A  simple  molecular-orbital  picture  of  the  wave  func¬ 
tions  of  the  £,£  states  and  of  the  H2  states  created  in 
channels  < la)  — ( lc)  suggests  that  each  of  these  channels 
can  be  related  to  specific  one-electron  transitions.  The 
transitions  can  be  easily  assigned  by  considering  the  form 
of  the  initial  and  final  electronic  states  for  each  channel. 
Near  the  E  well,  the  E,F  electronic  wave  function  is 
predominantly  \sog2scrg,  and  direct  ionization  of  the  E 
state  corresponds  to  ionization  of  the  2 sog  orbital. 
Channel  (lb1  corresponds  to  core  excitation  of  the  lso-g 
orbital  to  create  the  2sog2pou  autoionizing  state. 
Near  the  F  well,  the  electronic  wave  function  is  a  linear 
combination  of  \s<jg2sog  and  (2p<7u  )2.  In  this  range,  the 
autoionizing  state  can  also  be  formed  by  the 
2 pou  -~2sot  transition. 

In  the  third  process  (lc)  atomic  fragments  are  formed 


following  one-photon  absorption  of  the  E,F  level.  Direct 
dissociative  photoionization  to  the  continuum  of  the 
ground  state  and  first  excited  state  of  H2  +  is  also  possible 
from  the  F  well.  The  electronic  wave  function's 
\sag2sou  component  permits  coupling  to  the  ~2g 
ground  state  and  the  (2pcr„)2  component  to  the  excited 
state  of  22* .  The  various  channels  shown  will  occur  with 
amplitudes  that  are  weighted  according  to  the  radial  dis¬ 
tribution  of  the  E,F  vibrational  wave  (unction. 


B.  Photoelectron  spectra 

We  now  discuss  our  results  for  specific  intermediate  vi¬ 
brational  states  of  the  £,£  potential.  We  use  the  notation 
vE  F  to  denote  the  vibrational  quantum  number  according 
to  the  Dressing  convention,  which  combines  £-  and  F- 
well  states.  To  indicate  states  localized  predominantly  in 
the  inner  or  the  outer  well,  we  use  vE  and  vF,  respective¬ 
ly.  The  inner  £-well  potential  curve  is  approximately 
parallel  to  the  H2*  ground  state  2lg.  Therefore  direct 
photoionization  of  £-well  states  is  expected  to  have  a 
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large  component  from  transitions  with  |p  +  —  u£|=0.  An 
example  for  this  is  the  vibrational  level  vE  F  —  3  for  which 
the  wave  function  is  almost  entirely  localized  in  the  E 
well  (ii£  =  1).  The  most  likely  final  H2+  state  in  this  case 
is  v  *  =  1,  as  shown  in  Fig.  2.  Similarly,  as  shown  in  Fig. 
3,  d££  =  6  (u£  =  2)  leads  to  a  large  proportion  of  ti  +  =2. 
However,  in  both  cases  the  population  of  the  other  vibra¬ 
tional  states  is  larger  than  one  would  expect  from  simple 
Franck-Condon  arguments.  This  result  agrees  well  with 
the  result  of  a  similar  experiment  by  Anderson,  Kubiak, 
and  Zare,1  which  stimulated  theoretical  activity  to  ex¬ 
plain  the  anomalous  final-state  distribution.  We  note  that 
Anderson,  Kubiak,  and  Zare1  did  not  observe  the  v  ”  =  13 
peak  in  their  electron  spectrum  on  the  vEF=t  state, 
which  may  have  been  caused  by  a  low  sensitivity  in 
measuring  less  energetic  photoelectrons  in  their  electron 
collector. 

Our  theoretical  results  for  vEF  =  l  and  6  are  also 
shown  in  Figs.  2  and  3.  The  calculations  allow  us  to  as¬ 
sess  how  much  of  the  total  photoionization  comes  from 
the  direct  process  and  how  much  from  the  indirect  pro- 
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FIG  2.  Experimental  photoelectron  spectrum  for  the 
if  r  =  i  trt  =1’.  7=1.  intermediate  state,  and  the  theoretical 
distribution  of  final  vibrational  states.  The  shaded  bars  corre¬ 
spond  to  the  calculated  direct  cross  sections,  and  the  open  bars 
to  the  sum  of  the  direct  and  indirect  cross  sections.  The  in¬ 
direct  contribution  corresponds  to  the  difference  between  the 
two  bars 
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FIG.  3.  Same  as  Fig.  2,  for  i’£  f  =  6  (v£  =  2). 


cess  through  the  autoionizing  state.  The  shaded  bars 
show  the  Franck-Condon  portion,  and  the  open  bars 
show  the  total  calculated  cross  section.  In  general,  the 
agreement  with  experiment  is  good,  and  confirms  the  fol¬ 
lowing  interpretation  of  the  data.  The  dominant  chan¬ 
nels  are  (la)  and  (lb),  with  the  excitation  to  the  autoioniz¬ 
ing  state  contributing  to  the  broad  distribution  of  final 
states.  However,  other  mechanisms  should  not  be 
neglected.  Rudolph  et  al.2]  have  shown  that  effects  such 
as  the  R  dependence  of  the  transition  matrix  element  also 
contribute  to  the  broadening.  There  is  also  a  significant 
component  of  direct  photoionization  of  i)£/-  =  6  leading 
to  e  +  =  13.  The  two  largest  Franck-Condon  factors  for 
these  E,F  states  are  to  u  +  =2  (0.754)  and  to  u  *  =  13 
(0.089).  The  latter  component  arises  because  the  small 
outer-well  component  of  the  d££  =  6  state  contributes  to 
a  substantial  Franck-Condon  factor  with  v  *  =  13  of  H-^. 
The  calculations  also  show  that  the  indirect  ionization 
channel  contributes  significantly  to  the  v  *  =  13  peak. 

The  predominantly  inner-well  vibrational  states  would 
not  be  expected  to  have  a  large  cross  section  for  dissocia¬ 
tive  ionization  channel  (lc)  because  of  the  poor  Franck- 
Condon  factor  at  the  wavelength  used  in  this  experiment. 
The  experimental  data  confirm  this  expectation.  In  a 
separate  apparatus  that  monitors  ion  species  using  the 
same  laser  intensity  and  frequency,  H~  was  observed 
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only  as  a  small  percentage  of  the  total  ion  production. 
This  conclusion  agrees  with  Anderson,  Kubiak,  and 
Zare’s1  observation  on  the  i’E:F  =  6  state.  However,  Buck, 
Parker,  and  Chandler25  observed  a  somewhat  larger  pro¬ 
portion  of  H4. 

We  also  investigated  states  associated  with  the  outer,  F 
well.  The  u£f  =  7  (i>£  =  4),  J=  1  state  was  an  interesting 
choice  because  it  is  only  252.25  cm'1  higher  than  its  E- 
well  neighbor,  the  T££  =  6  (t>£  =  2),  J=  1  state.  The 
final-state  distributions  from  the  two  states  are  complete¬ 
ly  different.  Figure  4  shows  the  measured  spectrum. 
Several  low  v  4  states  have  comparable  signals,  and  the 
large  v  4  states  that  have  a  good  overlap  with  the  outer- 
well  wave  function  become  the  dominant  feature  in  the 
spectrum.  The  spectrum  from  vE  F  =  9,  shown  in  Fig.  5, 
also  exhibits  features  at  large  v  4  that  can  be  attributed  to 
the  overlap  between  the  outer-well  intermediate  state  and 
the  H:~.  This  electron  energy  spectrum  is  composed  of 
two  groups  of  peaks.  One  corresponds  to  the  state  with 
the  largest  Franck-Condon  overlap  with  the  inner-well 
portion  of  the  E,F  wave  function.  The  other  group  is 
dominated  by  v  4  =  14,  a  high  vibrational  state  from  ex¬ 
citing  the  outer  well.  We  also  notice  in  this  spectrum  a 
broad  electron  signal  below  1.5  eV  after  the  general  back¬ 
ground  was  subtracted  from  the  total  signal. 

Theoretical  calculations  for  t>££  =  7  and  9  are  also 
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FIG.  4  Same  as  Fig.  2.  for  vE  F  =  l  (vE  =4). 
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FIG.  5.  Same  as  Fig.  2,  for  vE  F  =  9  (r£  =  3). 


shown  in  Figs.  4  and  5.  The  interpretation  of  the  data  is 
not  as  clear  for  these  states  as  for  u££  =  3  and  6.  Because 
the  E,F  wave  functions  are  spread  over  both  wells,  the 
Franck-Condon  factors  for  the  direct  process  cover  many 
values  of  u  4 .  The  theoretical  calculations  indicate  that 
most  peaks  have  substantial  contributions  from  both 
direct  and  indirect  ionization  mechanisms. 

The  ion  mass  spectrometry  measurements  show  that 
ionization  via  the  vEF  =  9  (u£  =  3)  state  produces  a 
significant  amount  of  H4.  Figure  6  shows  the  TOF  spec¬ 
trum  of  ion  species  of  this  state;  the  H  4  signal  is  about 
30 7c  that  of  H24.  The  origin  of  this  large  fraction  of 
protons  will  be  discussed  in  Sec.  IV  C. 

C.  Origin  of  H  4 

Although  protons  account  for  a  significant  fraction  of 
the  photoions  observed  for  vE  F  =  9  in  the  present  experi¬ 
ment  and  in  other  recent  studies,  no  definitive  explana¬ 
tion  of  their  origin  exists.  Several  alternative  mecha¬ 
nisms  are  possible.  Referring  back  to  the  available  chan¬ 
nels  (la)-(lc),  we  note  that  excited  hydrogen  atoms  H* 
or  molecular  ions  H:4  can  be  produced.  Photoionization 
of  H*  or  photodissociation  of  H24  can  produce  protons. 
Channel  (lc)  produces  protons  directly  by  dissociative 
ionization.  Both  the  224  and  the  224  potential  curve  of 
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FIG.  7.  Results  of  rate-equation  model  for  the  fraction  of 
H2’  ions  that  are  photodissociated  to  form  H". 


FIG.  6.  Time-of-flight  spectrum  for  the  ion  mass  spec¬ 
trometry  measurements. 


H2+  are  available  for  this  process.  The  present  photo¬ 
electron  spectra  provide  additional  information  to  ad¬ 
dress  this  point. 

The  data  allow  us  to  rule  out  photoionization  of  H*  us¬ 
ing  a  fourth  photon  from  the  laser  pulse  as  a  significant 
source  for  the  observed  protons.  Excited  atoms  are  pro¬ 
duced  by  those  molecules  in  the  autoionizing  state 
that  "survive”  the  autoionization  process  and  dissociate. 
When  the  excited  atoms  produced  by  the  dissociation  are 
ionized,  the  photoelectron  energies  are  3.129  eV  for  n  —  2, 
5.017  eV  for  n  =  3,  and  5.678  eV  for  n  = 4  (at  the  laser 
wavelength  employed  for  vE%F=9).  We  have  not  ob¬ 
served  electrons  of  these  energies.  If  they  were  present, 
they  would  be  observed,  as  they  were  in  our  previous  ex¬ 
periment6  involving  photoionization  of  H2  (C  'll*  ). 

We  have  also  estimated  the  production  of  protons  by 
photodissociation  of  H2+,  using  a  rate-equation  model. 
Let  N0,  NE  F,  N  *,  and  denote  the  number  densities 
of  H2(A"),  H2(£,£),  H2*(u  * ),  and  protons,  respectively. 
Let  <7p,(i’*)  be  the  cross  section  for  photoionizing  the 
£,£state  and  producing  H2*(u  + ),  and  let  <7pd(v  + )  be  the 
cross  section  for  photodissociating  H2+(u  + ).  api(v  + )  de¬ 
pends  on  the  vibrational  state  of  the  £,£  potential  and 
Opjtu*)  depends  on  wavelength.  We  define  also 
<rpl  =  5;crp,(i>  *  1,  the  total  photoionization  cross  section  of 
the  £,£  state  and  JVH  -,  the  total  H2+  population.  The 
two-photon  excitation  cross  section  is  o0  and  has  units  of 
(length  i4  X  (time).  We  write  the  following  rate  equations, 
assuming  only  that  the  laser  pulse  is  not  strong  enough  to 
deplete  the  ground-state  population  A'0: 

dS  E  F  ^  ^  (  ^ 

—  =i'a0.\0  —  d>ar,.’^E  F  , 

d\H- 


These  equations  involve  explicitly  only  the  total  H2+ 
ion  population.  The  dependence  on  the  specific  popula¬ 
tion  of  a  vibrational  level  u  +,  is  contained  in 
which  is  a  weighted  average  over  the  photodissociation 
cross  sections  for  the  individual  vibrational  levels: 

^==2^  +  •  (3) 
2 

A  useful  estimate  of  the  proton  population  can  be  ob¬ 
tained  by  solving  the  above  equations  analytically,  using 
a  constant  value  for  a^.  We  can  obtain  an  estimate  of 
Opj  by  noting  that  the  H2+  states  are  populated,  at  least 
initially,  in  proportion  to  the  cross  sections  <7pj.  We 
therefore  assume  Nv  ♦  /NH  +  *apj(u  +  )/api. 

The  set  of  equations  (2)  with  constant  coefficients  is 
easily  solved  using  Laplace  transforms.  The  important 
dimensionless  parameters  are  (i)  ^opjt,  where  <t>  is  the 
laser  fluence  in  photons  per  cm_2s_l,  and  (ii)  0^/6^. 
Figure  7  shows  the  calculated  results  for  the  fraction  of 
total  ions  (H2+  and  H+)  that  are  protons.  The  rate  of 
converting  the  H2+  to  protons  by  photodissociation  com¬ 
petes  with  the  production  of  new  H2+  by  photoioniza¬ 
tion;  the  balance  depends  on  the  ratio  a^/a^. 

Typical  values  of  cr^/a^  are  in  the  range  0.4-0.6. 
These  values  were  determined  from  the  present  calcula¬ 
tions  of  op)  and  from  direct  calculations  of  the  photodis¬ 
sociation  cross  section  similar  to  Dunn’s.24  From  Fig.  7 
it  is  clear  that  for  photodissociation  of  H2*  to  contribute 
substantially  (  ~  10%)  to  the  proton  fraction,  <bopit  —  1  is 
necessary.  Typical  parameters  for  the  present  experiment 
are  100  /iJ  per  5  ns  pulse,  and  spot  size  500  /im.  These 
values,  combined  with  the  calculated  values  of 
crp,~7X  10" 18  cm:,  yield  ar  estimate  of  <bapti  ~0.1. 

The  rate-equation  analysis  leads  to  the  conclusion  that 
photodissociation  of  is  likely  to  account  for,  at  most, 
a  few  percent  of  the  observed  protons  in  the  present  ex¬ 
periment.  We  note  that  substantially  more  protons 
would  be  expected  from  this  mechanism  at  greater  laser 
intensities. 

Another  mechanism  for  producing  protons  is  dissocta- 
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tive  ionization  of  the  E,F  states.  Continuum  translation¬ 
al  states  on  both  the  2I  +  and  21^H2  +  potential  curves 
are  possible  final  states  for  this  process.  Dissociative  ion¬ 
ization  can  also  occur  following  excitation  of  the  '2*  au- 
toionizing  state.  Here  we  discuss  qualitatively  the  direct 
photoionization  to  continuum  states  of  the  21+  and  22* 
potentials.  The  cross  section  for  this  process  can  be  writ¬ 
ten25 

=  ~1  |J  Y/(K  >M(*  >*,<*><**  f  •  W) 

where  is  the  initial  vibrational  wave  function  on 

the  E,F  curve,  fj.(R)  is  the  transition  dipole-matrix  ele¬ 
ment,  and  )  is  the  translational  wave  function  on  the 
H  +  H  +  potential  curve.  The  energy  supplied  by  the  pho¬ 
ton  is  shared  by  the  photoelectron  and  the  nuclei;  one 
must  calculate  the  photoelectron  spectra  expected  by 
evaluating  Eq.  (4)  for  wave  functions  >  for  *he  range 
of  H  +  H+  energies  possible.  Calculations  of  fiiR)  are 
not  available,  but  qualitative  conclusions  follow  from 
comparing  the  Franck-Condon  factors  for  transitions  to 
the  final  states  21*  and  22*.  Direct  numerical  integra¬ 
tion  of  Eq.  (4)  for  fi=  I  and  for  photon  energies  corre¬ 
sponding  to  Vg'F  =  9  shows  that  the  Franck-Condon  fac¬ 
tors  to  the  continuum  of  the  22~  state  are  large  only  for 
nuclear  translational  energy  states  in  the  energy  range 
S 100  meV.  For  large  translational  energies,  the 
Franck-Condon  factors  greatly  favor  the  2i*  final  state. 

An  estimate  of  the  absolute  contributions  of  dissocia¬ 
tive  ionization  via  the  continuum  of  the  ground  state  of 
H2*,X2Ig  can  be  obtained  by  summing  the  calculated 
contributions  of  direct  photoionization  to  all  bound  levels 
of  H2*.Y.  In  Table  I  we  give  these  sums  for  the  vibration¬ 
al  levels  investigated  here.  It  is  apparent  that  only  for 
vE  F  =  9  does  a  substantial  contribution  from  the  ground- 
state  continuum  exist  i  —69c).  It  is  most  likely  that  this 
contribution  together  with  the  transition  to  22*  are  re¬ 
sponsible  for  the  proton  yield  observed  in  vE  r  =  9. 

The  processes  tic  must  give  rise  to  photoelectrons  at 
energies  below  1.53  eV  in  Fig.  5.  The  data  suggest  that 
electrons  are  detected  which  correspond  to  this  dtssocia- 


TABLE  1.  Summed  Franck-Condon  factors  to  the  ground 
state  of  H;*  (calculated  for  7=0). _ 


»e.f 

2  l<»£>*>lJ 

r  +  -0 

3 

1.000 

6 

0,999 

7 

0.989 

9 

0.942 

tive  ionization  process,  although  more  refined  experi¬ 
ments  will  be  required  to  prove  the  undulatory  electron 
energy  distributions  predicted  by  Eq.  (4). 

V.  CONCLUDING  REMARKS 

We  have  reported  measurements  and  calculations  of 
the  photoelectron  energy  spectra  following  resonant  mul¬ 
tiphoton  ionization  of  H2  involving  intermediate  vibra¬ 
tional  states  on  the  double-well  E,F  potential  curve.  The 
available  ionization  channels  include  direct  photoioniza¬ 
tion,  indirect  photoionization  through  a  dissociative  au- 
toionizing  state,  and  dissociative  ionization.  Our  results 
show  that  the  relative  importance  of  these  channels  de¬ 
pends  sensitively  on  whether  the  intermediate  E,F  wave 
function  is  localized  in  the  E  or  the  F  well.  In  the  latter 
case,  qualitative  arguments  based  on  Franck-Condon  fac¬ 
tors  support  the  interpretation  that  most  of  the  dissocia¬ 
tive  ionization  occurs  through  fhe  22^  state.  Another 
likely  source  of  H+  production,  photodissociation  of 
H2*.  does  not  appear  to  be  a  dominant  mechanism  at  the 
present  power  levels. 
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ABSTRACT 


Calculations  arc  reported  of  the  dissociative  attachment  of  low  energy  (<  5eV)  electrons  to 
molecular  hydrogen  in  the  states  v  =  0-9  and  J  *  0.  The  dynamics  arc  treated  using  resonant 
scattering  theory,  fully  including  the  nonlocal  and  energy-dependent  shift  and  width  operators. 

The  necessary  coupling  terms  are  taken  from  the  ab  initio  calculations  of  fixed-R  electron-^ 
scattering  performed  by  Miindel,  Berman,  and  Domcke.  Good  agreement  with  the  available 
experimental  data  of  Schulz  and  Asir>  li  and  of  Allan  and  Wong  is  achieved.  Exploratory 
calculations  were  performed  to  investigate  various  local  approximations  to  the  nonlocal  shift 
operator.  The  results  suggest  that  the  optimum  local  approximation  will,  in  general,  differ  from  the 
curve  obtained  in  a  fixed-nuclei  structure  calculation. 
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INTRODUCTION 


Dissociative  attachment  (DA)  of  electrons  to  molecular  hydrogen  is  a  process  of  great 
fundamental  and  practical  interest.  It  is  thought  to  be  a  primary  source  of  the  H~  ions  produced  in 
hydrogen  plasmas.  The  process  may  be  schematically  written 

e  +  H2(v,J)  H’+  H  (1) 

where  v  and  J  denote  the  vibrational  and  rotational  quantum  numbers  of  the  H2.  Measurements  for 
the  absolute  cross  section  for  attachment  to  the  ground  state  (v  =  J  =  0)  were  reported  by  Schulz 
and  Asundi,1  and  Allan  and  Wong2  showed  that  the  cross  section  increases  dramatically  as  the 
internal  energy  of  the  H2  increases.  They  measured  relative  values  of  sigma  for  v  =  0-4  (J  =  0), 
and  found  that  the  cross  section  increased  roughly  an  order  of  magnitude  for  each  increment  of  v. 

This  behavior  was  explained  in  1978  by  the  analysis  of  Wadehra  and  Bardsley3  based  on 
the  local  complex  potential  model.  They  used  ab  initio  potential  curves  available  at  that  time,  and 
fit  a  parameterized  coupling  function  to  the  experimental  data.  The  same  coupling  function  was 
then  used3^  to  extend  the  calculations  to  other  values  of  (v  J),  for  which  experimental  data  are  not 
available. 

Recent  work  in  the  calculation  of  potential  curves  and  coupling  terms  for  electron  molecule 
scattering  has  provided  ab  initio  information  of  much  greater  reliability  than  was  available  ten  years 
ago,  and  this  information  can  be  used  to  reexamine  the  treatment  of  dissociative  attachment.  For 
example,  the  recent  calculations  of  Mundel,  Berman  and  Domcke7  have  provided  the  necessary 
matrix  elements  for  a  treatment  of  dissociative  attachment  of  H2  that  fully  includes  the  energy 
dependent  and  nonlocal  shift  and  width  operators.  These  authors  also  report  calculations8  of  DA 
of  H2  for  v  =  0-2  (J  =  0),  but  the  emphasis  of  their  work  was  on  vibrational  excitation.  We  have 
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found  that  their  work,  with  some  additional  analysis,  provides  a  basis  for  much  broader 
calculations  of  DA. 

This  paper  is  organized  as  follows:  Section  II  presents  the  theory,  including  a  summary  of 
the  resonant  scattering  theory  and  a  discussion  of  specific  coupling  terms.  Particular  care  is  taken 
to  handle  the  potentials  and  couplings  in  the  range  of  intemuclear  separations  R  near  3.0  a<,  where 
the  resonant  state  H2  crosses  the  H2  potential  curve.  Section  HI  presents  and  discusses  the 
numerical  results,  and  Section  IV  contains  concluding  remarks. 
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THEORY 


A.  SUMMARY  OF  RESONANT  SCATTERING  THEORY 

The  formal  theory  of  DA  is  well  developed,8"11  and  we  will  present  only  the  essential 
points  to  facilitate  discussion  of  our  implementation.  The  following  summary  follows  the  notation 
of  Wadehra's  review12  to  a  large  extent.  We  begin  with  the  expression  for  the  total  wave  function 
of  the  electron-molecule  system: 

'F(q,R)  =  <t>d(q,R)F(R)+  I  J  de  fv(e)  ye  <q,R)  Xv  (R)  (2) 

v 

open 

q  refers  collectively  to  the  electron  coordinates,  and  R  is  the  intemuclear  separation.  It  is  assumed 
that  the  fixed-R  electron-molecule  wave  functions  have  been  determined:  <j>d  is  a  discrete  state, 
corresponding  to  Hjj  and  \jre  is  an  electron-molecule  scattering  wave  function.  The  vibrational 
wave  functions  of  H2  are  Xv-  The  treatment  of  DA  requires  the  determination  of  F(R),  the  wave 
function  describing  the  nuclear  motion.  By  substituting  Eq.(2)  into  Schroedinger’s  equation  and 
formally  solving  for  the  expansion  coefficients  fv(R),  one  can  derive  the  following  equation  for 
F(R): 

[Iff  -  Vd(R)  +  e]  F(R)  -  V(R,E-EVo)  XVo(R)  +  (A  -  jO  F(R)  (3) 

The  total  energy  E  is  the  sum  of  the  asymptotic  (R  -» <»)  kinetic  energy  of  the  electron  and  the 
vibrational  energy  of  the  molecule,  which  is  Ev  for  vibrational  state  v.  The  initial  vibrational  state 
of  H2  is  Vq.  V<j(R)  is  defined  as  a  matrix  element  of  the  electronic  Hamiltonian  He] 

Vd(R)  =  «bd 1  Hei  I  0d>  (4) 
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A  and  T  are  nonlocal,  energy-dependent  operators  whose  action  on  the  function  F(R)  yields  a 
function  of  R  defined  as  follows: 


V(R,  E)  =  «t>d  I  Hei  I  Ve> 


(8) 


Equation  (3)  is  solved  for  the  boundary  conditions  that  F(R)  =  0  for  R  =  0,  and  F(R)  is  a 
outgoing  wave  asymptotically.  The  cross  section  for  DA  is  then 


ODA--1J  Bm  IF(R)P 

e  M  R— 


This  formulation  can  be  related  to  the  parameters  describing  the  fixed-R  resonant  scattering. 
At  each  R,  the  electron  scattering  cross  section  from  the  molecule  exhibits  a  resonance  of  width 
F(R)  at  an  electron  kinetic  energy  Eres(R)-  It  is  convenient  to  define  a  potential  curve  V^OR)  that 
combines  the  electron  scattering  energy  with  the  energy  of  the  bound  state  H2  target: 

Vres(R)  =  VH2(R)  +  Eres(R)  (10) 
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where  Vh2(R)  is  the  ground  state  H2  potential.  Vres(R)  is  shifted  by  an  amount  A(R)  from  the 
zeroth  older  discrete  state  energy  Vd(R).  Y^fR),  E^fR)  and  T(R)  are  related  as  follows: 


Vres(R)  =  Vd(R)  +  A(R)  (11) 

where 

A(R)  =  A(R,R,Eres(R))  (12) 

The  width  (R)  of  the  resonance  at  fixed  R  is  given  by 

F(R)  =V(R,Eres(R))2  (13) 

The  distinction  between  V^R)  and  Vd(R)  has  not  always  been  maintained.  It  is  often 
assumed  that  the  operator  A  on  the  right  hand  side  (rhs)  of  Eq.(3)  can  be  neglected,  and  that  Vd(R) 
on  the  lhs  can  be  replaced  by  V^CR) .  This  simplification  corresponds  to  approximating  the 
dynamical  operator  A  by  A(R)  in  Eq.(3),  and  then  assuming  the  closure  relation  for  the  sum  over 
open  vibrational  states.  In  the  present  work,  we  utilize  recent  calculations7*8  for  V(R,e),  which 
allows  us  explicitly  to  include  the  nonlocal  shift  operator. 

B  .  ELECTRONIC  MATRIX  ELEMENTS 

Calculations  are  available  in  the  literature  for  nearly  all  the  necessary  coupling  terms  for  DA 
of  H2.  Excellent  potential  curves  for  H2  have  been  obtained13;  these  were  used  to  calculate  the 
vibrational  wave  functions  Xv.  E^R)  has  been  calculated  by  Mundel  et  al.7  These  authors 

showed  that  Eres(R)  could  be  unambiguously  defined  as  the  electron  scattering  energy  at  which  the 
eigenphase  sum  equals  tt/2. 

For  the  essential  matrix  element  V(R,e),  we  have  fit  the  numerical  values  calculated  by 
Miindel  et  al.7-8  by  a  convenient  analytic  function.  Although  Mtlndel  et  al.8  reported  a  nine- 
parameter  analytic  fit  to  their  calculations,  we  have  found  it  preferable  to  repeat  the  fit  using  a 
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different  analytic  form.  Some  difficulty  was  encountered  in  trying  to  use  the  functional  form 
reported  in  Ref.  8.  V(R,e)  approaches  a  nonzero  constant  asymptotically  (as  R— *»).  Also, 
Table  I  in  Ref.  8  contains  typographical  errors. 

A  good  fit  to  the  calculated  points  was  obtained  with  the  following  analytic  function: 


V(R,e)  = 


A(R)e3# 
(e  +  B)2 


(14) 


This  formula  is  physically  plausible.  If  we  consider  the  matrix  element  between  a  bound 
hydrogenic  function  and  an  energy-normalized  p  wave  scattering  function  ^k/rc)1/2  ji(kr),  the 
result  is14 

oo 

<2k/n)>'2  I  e-"'°  ji  (kr)  r2  dr  ~  (~2  *”.2)2  (15) 

Since  e  =  (l/2)k2,  Eq.(15)  is  consistent  with  Eq.(14)  and  yields  the  correct  threshold  behavior. 

The  most  accurate  fit  could  be  achieved  by  letting  both  A  and  B  be  functions  of  R  in 
Eq.(14).  However,  considerable  simplification  of  the  calculations  results  if  B  is  a  constant.  This 
point  will  be  discussed  below.  Fortunately,  we  found  that  a  reasonably  accurate  fit  could  still  be 
obtained  if  B  is  a  constant.  Fig.  1  shows  the  results  of  a  least  squares  fit  to  the  calculations  of 
Miindel  et  al.,7  using  the  analytic  form  of  Eq.(14).  The  constant  B  =  7.15  eV  was  obtained  in  the 
fit  This  approach  also  provides  a  natural  way  to  extrapolate  V(R,e)  to  large  values  of  R:  one 
keeps  the  same  B  and  adopts  a  sensible  asymptotic  form  for  A(R).  As  shown  in  Fig.  2,  the  values 
of  A(R)  fit  in  the  range  1.4  ao  <  R  <  2.75  ao  (for  which  calculations  were  reported)  suggest  a  linear 
extrapolation  out  to  the  R  for  which  A(R)  is  zero. 

The  advantage  of  the  form  Eq.(14)  for  V(R,e)  is  twofold.  First,  the  principal  part  integral 
in  the  definition  of  the  shift  operator  [Eq.(6)]  can  be  done  analytically.  Second,  a  separable  form 
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of  the  shift  operator  results,  permitting  Eq.(3)  to  be  solved  by  Bardsley’s  Green  function 
method.15-16  We  consider  the  principal  part  integral  in  this  section,  and  defer  discussion  of  the 
method  of  solution  to  the  next  section. 

We  now  consider 


A(R,R',e)  =  P 

0J  (e'+B)4  (£-£’) 


(16) 


where  P  denotes  the  principal  part  integral.  Simplification  follows  by  change  of  variable  to  x  = 
£'/B,  and  by  the  introduction  of  a  factor  of  ti/16  for  future  convenience: 


A(R,R',£)  =  A(R)  A(R')  B5/2  P(e/B) 


(17) 


where 


7t 


OO 

Pf 

o 


X3/2  dx 
(x+1)4  (x-a) 


(18) 


The  function  p(a)  only  needs  to  be  evaluated  in  the  interval  [-1,1].  The  value  of  p  elsewhere 
follows  from  the  symmetry  relation 


p(l/a)  =  -a  p(a),  (19) 

which  is  easily  verified  by  writing  the  integral  expression  Eq.(18)  for  p(l/a),  and  then  changing 
variables  to  y  =  1/x. 
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Fig.  3  illustrates  the  values  of  p(a)  calculated  numerically.  The  special  values  p(0)  =  1, 


p(l)  =  0,  and  dp/da  (0)  =  5  follow  from  Eqs.(18)  and  (19). 

We  now  illustrate  the  explicit  values  of  the  coupling  matrix  elements  used.  Particular 
attention  is  paid  to  the  region  of  R  near  3  a<>  For  R  >  3  ao,  there  exists  a  well  defined  bound  state 
Hi  at  an  energy  below  the  corresponding  H2  energy.  For  fixed  R  <  3  ao,  there  exists  a  scattering 
resonance  at  an  energy  Eres(R)  >  0.  Formally,  one  expects17  the  bound  state  at  large  R  to  join 
smoothly  with  the  resonance  at  small  R;  the  potential  curves  adopted  to  treat  DA  must  behave  this 
way  as  well.  Using  the  values  of  Eres(R)  calculated  bu  Mundel  et  al.7  for  R  <  3  ao,  and  the 
potential  curve  calculated  by  Senekowitsch  et  al.18  for  R  >  3  ao,  we  found  indeed  that  the  Hi  curve 
at  large  R  could  be  smoothly  joined  with  the  curve  Vres(R)  defined  by  Eq.(10).  Another  way  of 
stating  this  is  that  for  R  >  3  ao,  Eres(R)  is  negative  and  corresponds  to  a  bound  state. 

We  can  now  define  the  potential  curve  Vd(R)  for  all  values  of  R.  Since  EreS(R)  has  been 
shown  to  be  a  smooth  and  well-behaved  function  of  R  near  R  =  3  a<j,  there  is  no  difficulty  in  using 
the  definition  of  the  shift  [Eqs.  (11),  (12),  and  (17)]  to  determine  Vd(R)  for  both  positive  and 
negative  values  of  Eres(R): 


Vd(R)  =  Vh2(R)  +  Eres(R)  -  A  (R)  (20) 

The  potential  curves  determined  are  shown  in  Fig.  4. 

C.  WORKING  EQUATIONS 

We  now  describe  the  solution  of  Eq.(3).  The  analysis  of  the  previous  section  has  led  to 
explicit  forms  of  the  shift  and  width  operators  that  can  be  handled  with  standard  methods.  The 
crucial  point  is  that  the  R  and  R'  dependence  of  A  and  T  is  identical.  One  can  rewrite  Eq.(3)  as 
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where 


and 


^2  -  Vd(R)  +  e]  F(R)  =  V(R,  E-EVo)  XVo  (R) 

+  J  dR'  I  Xv  (R)  Xv  (R’)  A(R)  A(R') 
v 

open 

x  [gA  (E-Ev)  -  in  gp  (E-Ev)]  F  (R1) 


(21) 


gr  (e)  = 


£3/2 

(e  +  B)4 


(22) 


gA  (e)  =  -^B-5/2p(e/B) 


(23) 


Bardsley's  Green  function  method15-16  can  be  used  to  reduce  the  problem  to  a  set  of  linear 
equations  whose  coefficients  involve  the  matrix  elements  of  the  Green  function  for  the  lhs  of 
Eq.(20).  This  function  is  defined  as 

G(R',R)  =  C’1  F,(R<)  F2(R>)  (24) 

where  Fi(R)  is  the  regular  solution  of  Eq.(20)  with  the  rhs  set  to  zero  that  satisfies  the  boundary 
conditions  Fi(0)  =  0  and  Fi(R)  — »  sin(kR+8)  as  R-*».  The  irregular  solution  F2(R)  satisfies 

F2(R)-»  exp[i(kr+8)]  as  R-»«>.  For  this  choice  of  asymptotic  normalization,  the  Wronskian 
C  =  -k/2M.  The  matrix  elements  G  .  are  defined  by 

oo  oo 

Gvv.  =  J  dR  J  dR'  A(R)Xv(R)  G(R,R,)A(R’)XV-(R')  (25) 

o  o 

One  then  solves  the  set  of  linear  equations  Ax  =  b,  where 
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oo 


(26) 


xv  =  /  XV(R)  A  (R)  F  (R)  dR 
o 

bv  =  gf  (E-EVo)  GWo  (27) 

and 

Ayv'  =  Syv'  •  £g^(E*Ev')  -  IK  gp(E-Ey')^ J Gyy'  (28) 


The  asymptotic  form  of  F(R)  can  be  determined  from  the  xv  just  as  in  previous  work,16-19  and 
from  that  the  cross  sections  can  be  evaluated. 
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RESULTS  AND  DISCUSSION 

A.  CALCULATIONS  OF  CROSS  SECTIONS  AND  RATES 

Cross  sections  for  DA  of  H2  through  the  22^  H2  state  have  been  calculated  for  J  =  0  and 
v  =  0-9,  and  for  electrons  from  threshold  up  to  4-5  eV.  Results  for  the  cross  sections  as  a 
function  of  incident  electron  energy  are  shown  in  Fig.  5.  These  calculations  include  the  full 
energy-dependent  and  complex  shift  and  width  operators.  The  cross  sections  exhibit  oscillatory 
structure  near  threshold  that  appears  to  be  related  to  the  energy  dependence  of  the  width  operator. 
We  found  that  the  oscillations  persisted  when  the  nonlocal  shift  was  replaced  by  a  local  shift. 
However,  the  structure  was  not  present  in  earlier  calculations  we  performed  with  an  nonlocal  but 
energy-independent  width,  nor  was  it  seen  in  the  calculations  of  Wadehra  and  Bardsley,3  who 
used  a  local  shift  and  a  local  width.  This  comparison  suggests  that  the  energy-dependent  factors  in 
the  sum  over  vibrational  states  in  Eq.(21)  are  the  origin  of  the  structure.  These  terms  spoil  the 
closure  approximation,  in  which  the  sum  is  replaced  by  a  Dirac  delta  function,  and  which  underlies 
the  reduction  to  the  local  complex  potential  model. 

A  more  physical  explanation  is  the  following:  The  local  width  approximation  means  that 
the  resonant  H2'  state  decays  at  a  smoothly  varying  rate  as  the  nuclei  move  apart  after  electron 
capture.  This  decay  rate  is  independent  of  the  details  of  the  final  channels  available.  The  nonlocal 
width  restricts  the  decay  of  the  resonant  state,  by  introducing  a  dependence  on  the  matrix  element 
of  that  state  with  the  wave  functions  of  the  available  vibrational  levels  of  H2.  The  decay  is  then 
sensitive  to  the  details  of  the  matrix  elements  between  the  vibrational  levels  of  H2  and  the 
translational  wave  function  on  the  H2'  state.  Because  the  matrix  elements  change  with  energy,  the 
oscillatory  structure  results.  This  structure  may  also  be  an  analog  of  the  vibrational  fine  structure 
discussed  by  Domcke  et  al.1 1  in  e-H2  vibrational  excitation. 
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The  calculations  can  be  compared  with  experimental  data,  and  with  other  calculations  in  the 
literature,  in  several  different  ways.  Fig.  6  shows  the  cross  section  for  attachment  by  H2  (v  =  0) 
as  a  function  of  energy.  Several  recent  calculations  are  compared  with  the  data  of  Schulz  and 
Asundi.1  The  calculated  peak  values  are  very  similar.  It  is  plausible  that  the  calculations  are  all 
larger  than  the  experiment;  Gauyacq20  has  pointed  out  that  convoluting  his  calculation  with  an 
electron  energy  distribution  whose  full  width  at  half  maximum  is  300  meV  would  decrease  the 
calculations  from  2.3  to  1.5  x  10-5  A2.  The  primary  difference  between  the  calculations  is  the 
width  of  the  threshold  peak.  The  present  calculation  agrees  closely  with  that  of  Gauyacq,20  and 
obtains  a  FWHM  of  about  340  meV;  the  calculated  widths  of  Wadehra  and  Bardsley3  and  Mundel 
et  al.8  are  about  500  and  530  meV,  respectively.  Schultz  and  Asundi1  reported  two  sets  of  data 
taken  under  different  conditions.  The  set  with  the  narrower  width,  about  400  meV,  is  shown;  the 
width  of  the  other  set  was  about  600  meV.  Further  experimental  data  would  be  desirable  to 
differentiate  between  the  various  calculations. 

Allan  and  Wong2  have  measured  the  relative  cross  sections  at  threshold  for  v  =  0-4,  and 
we  compare  our  calculations  with  their  data  in  Fig.  7.  The  relative  experimental  measurements 
have  been  converted  to  absolute  cross  sections  by  using  the  measurements  of  Schulz  and  Asundi 
for  v  =  0  (=  1.6  x  10'5  A2).  Comparision  with  other  calculations  is  also  possible  over  the  range 
v  =  0-9.  We  note  that  Wadehra  and  Bardsley3'4  and  Wadehra5>6  fit  their  calculations  to  the 
experimental  data,  and  used  the  same  coupling  potentials  and  matrix  elements  for  the  calculations  at 
higher  v.  The  present  calculations  are  closest  to  those  of  Mundel  et  al.8  for  v  =  0-2,  who  also 
took  into  account  the  nonlocal,  energy-dependent  shift  and  width  operators.  The  present 
calculation  rests  on  a  fully  numerical  solution  of  the  nuclear  wave  equation  [Eq.(21)],  which  was 
facilitated  by  a  particular  choice  of  the  form  of  the  electronic  coupling  matrix  element  [Eq.(14)]. 
Mundel  et  al.  fit  Vd(R)  with  a  Morse  potential,  and  treated  the  shift  and  width  operators  by  an 
expansion  in  a  Lanczos  basis.  We  believe  that  the  present  approach  uses  more  accurate  potentials 
for  the  range  of  intemuclear  distance  where  the  H2  potential  crosses  the  H2  curve.  This  range  is 


crucial  to  DA,  but  less  important  to  vibrational  excitation,  which  was  the  emphasis  of  the  work  of 
Miindel  et  al.8  Finally,  we  note  that  the  present  calculations  agree  well  with  those  of  Gauyacq,20 
which  were  based  on  the  zero-range-potential  formalism. 

Extensive  calculations  were  carried  out  by  Wadehra  and  Bardsley3*4  and  Wadehra5*6  using 
the  local  complex  potential  model;  many  of  the  results  were  presented  as  rate  constants  depending 
on  an  average  energy  E  =  (3/2)kT.  We  adopt  the  same  nomenclature  for  easy  comparision  with 
their  results.  Fig.  8  shows  rather  close  agreement  between  the  present  calculations  and  the 
previous  calculations  based  on  the  local  complex  potential  model.  Although  the  logarithmic  scale 
minimizes  the  quantitative  differences,  all  general  trends  are  in  essential  agreement. 

Finally,  the  dynamical  calculations  of  DA  depend  on  the  potentials  and  matrix  elements  in 
the  range  of  R  >  3  aQ,  where  Hi  is  a  true  bound  state.  A  series  of  calculations  was  performed  in 
which  A(R)  was  arbitrarily  set  to  zero  for  R  >  3  ao,  thereby  eliminating  any  contribution  from  this 
range  of  R.  The  total  cross  sections  for  DA  decreased  about  20-30%.  Since  Miindel  et  al.7*8  only 
reported  numerical  values  for  V(R,e)  in  the  range  R  =  1.4  to  2.75  ao,  the  A(R)  we  used  for  larger 
R  is  an  extrapolation.  This  test  provides  a  measure  for  the  uncertainty  in  our  calculations  due  to 
our  extrapolation  of  A(R). 

B.  ACCURACY  OF  THE  LOCAL  COMPLEX  POTENTIAL  MODEL 

The  most  surprising  aspect  of  the  present  results  is  that  the  calculations  based  on  the  local 
complex  potential  (LCP)  appear  to  be  fairly  satisfactory,  at  least  when  one  looks  at  an  averaged 
quantity  like  the  rate  constant.  Is  there  a  physical  basis  for  replacing  the  nonlocal  operators 
A  -  (i/2)  T  with  a  simpler,  local  form?  Miindel  et  al  8  performed  a  comparison  for  v  =  0  and 
concluded  that  the  LCP  was  a  "stunning  failure".  They  noted  that  the  results  of  the  LCP  could  be 
as  much  as  a  factor  of  ten  large.  The  explicit  expressions  we  have  used  for  the  width  and 
particularly  the  shift  operator  provide  a  tool  for  monitoring  the  accuracy  of  the  final  calculation  as 
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successive  approximations  are  invoked  to  simplify  the  shift  or  width.  Our  results  confirm  the  need 
to  consider  the  full  operators,  but  suggest  situations  where  an  optimum  local  shift  can  be  defined. 

We  begin  our  analysis  with  the  explicit  expression  for  the  action  of  the  shift  operator  on  the 
nuclear  wave  function,  which  follows  from  Eqs.(21)  and  (23) 


AF  =  ^ B-5/2  A(R)  J  dR'  I  Xv  (R)  Xv  (R’>  p  A(R’)  F(R’)  (29) 

open 

The  operator  A  is  both  nonlocal,  because  of  the  integrals  over  R',  and  energy  dependent,  because 
of  the  terms  involving  (E-Ev).  We  note  that  the  constant  B  was  determined  to  be  7.15  eV,  and  that 
for  electrons  at  threshold  for  DA,  (E-Ev)  ranges  from  3.7  eV  (v  =  0)  to  0.14  eV  (v  =  9).  Hence  for 
our  calculations  (E-Ev)/B  <  Q.5.  From  Fig.  3  it  appears  that  the  energy  variation  of  p((E-Ev)/B)  in 
this  range  is  not  large  and  that  p(0)  =  1  could  be  substituted  into  Eq.(29)  as  an  approximation.  If 
this  is  done,  the  shift  operator  is  reduced  to  an  energy- independent  but  still  nonlocal  operator.  The 
arguments  for  further  simplification  are  well  known;  one  assumes  that  the  sum  over  open  channels 
in  Eq.(29)  is  close  to  the  sum  over  all  channels  (including  the  continuum),  which  satisfies  the 
completeness  relation 


I  Xv  (R)  Xv  (R')  =  S(R-R')  (30) 

v 

Again  specializing  to  H2,  we  note  that  the  sum  always  includes  at  least  10  terms,  since  the 
asymptotic  energy  of  H  +  H'  is  greater  than  the  vibrational  energy  of  H2  for  v  =  0-9.  If  we 
assume  that  the  sum  is  complete,  the  shift  operator  is  reduced  to  the  following  local  operator: 

AF  =  --j^  B-5/2  A(R)2F(R)  (31) 
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This  form  provides  a  simple  analytic  form  for  the  shift.  Determination  of  an  explicit  analytic  form 
for  the  coupling  matrix  element  V(R,e)  has  led  to  an  approximate  form  of  a  local  shift  that  is 
different  from  the  normal  local  approximation  to  the  shift  It  is  usually  argued  that  the  local  shift  is 
the  static  shift  that  is  obtained  when  one  solves  the  fixed  nuclei  scattering  problem  [Eq.(l  1)].  We 
have  tested  these  two  forms  of  the  local  shift. 

A  similar  argument  can  be  used  to  reduce  the  nonlocal  width  operator  to  a  local  form.  In 
this  case,  however,  there  is  no  clear  way  to  eliminate  the  energy  dependence  as  there  was  for  the 
shift  Substituting  Eres(R)  for  E-Ev  in  the  expression  for  r(R,R,E-Ev)  in  Eq.(5)  is  the  most 
straightforward  approach.  As  expected,12-21  test  calculations  with  such  an  energy-independent  but 
nonlocal  width  show  spurious  threshold  behavior:  the  cross  section  decreases  discontinuously 
(generally  less  than  10%)  as  a  new  channel  for  vibrational  excitation  opens  up.  The  calculations 
also  suggested  that  the  sum  over  open  vibrational  states  was  converged  within  10-20%. 

We  performed  a  series  of  calculations  to  test  the  various  approximations  to  the  shift 
operator.  In  all  cases  the  nonlocal  and  energy-dependent  width  was  preserved.  We  found  that  the 
local  shift  defined  by  Eq.(3i  was  a  reasonably  good  approximation.  Fig.  9  compares  the  results 
with  the  local  shift  and  with  the  full  shift  operator.  The  largest  error,  about  a  factor  of  two,  occurs 
for  v  =  0,  and  the  approximation  rapidly  converges  for  increasing  v.  This  behavior  is  completely 
consistent  with  the  arguments  that  led  to  the  local  approximation  in  the  first  place.  (E-Ev)/B  at 
threshold  is  largest  for  v  =  0,  and  becomes  progressively  smaller  for  larger  v.  Hence 
approximating  p((E-Ev)/B)  by  p(0)  should  become  progressively  more  accurate  for  larger  initial  v, 
and  this  behavior  is  exhibited  by  Fig.9. 

We  also  performed  calculations  in  which  the  local  shift  was  approximated  by  the  fixed-R 
shift  defined  by  Eq.(l  1).  As  before,  the  energy-dependent  and  nonlocal  width  operator  was  used. 
For  v  =  0  the  cross  sections  were  about  an  order  of  magnitude  too  large,  confirming  the 
observation  of  Miindel  et  al.8  that  this  form  of  the  local  shift  is  a  poor  approximation.  The  cross 
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sections  were  generally  less  accurate  than  those  obtained  with  Eq.(31),  although  the  results  were 
within  a  few  percent  for  v  =  5-9. 

The  two  local  potentials  corresponding  to  the  alternative  local  shifts  are  compared  in 
Fig.  10.  Viocai  is  obtained  by  adding  the  local  shift  defined  by  Eq.  (31)  to  Vd;  it  led  to  reasonably 
accurate  cross  sections.  Vres  was  defined  earlier  in  Eqs.  (11)  and  (12)  as  the  position  of  the  fixed- 
nuclei  electron  scattering  resonance.  It  is  the  curve  that  would  be  obtained  in  a  standard  fixed- 
nuclei  calculation  and  we  have  seen  that  it  does  not  lead  to  accurate  cross  sections. 

These  results  demonstrate  that  part  of  the  failure  of  the  LCP  reported  in  Ref.  8  can  be 
attributed  specifically  to  the  choice  of  a  local  shift,  and  that  the  LCP  can  be  improved  by  alternative 
choices  of  A(R).  However,  the  optimum  choice  in  the  present  case  is  not  the  fixed-nuclei  form  one 
might  expect  on  physical  grounds.  These  results  provide  a  rationale  for  the  success  of  Wadehra’s 
and  Bardsley's  calculations^-6  based  on  the  LCP  model.  Because  of  the  accuracy  of  a  local  shift 
constructed  by  well  defined  mathematical  approximations,  it  appears  plausible  that  a  fitting 
procedure  would  lead  to  a  legitimate,  approximate  local  shift  and  width.  However,  using  ab  initio 
fixed-nuclei  shifts  and  widths  in  an  LCP  model  calculation  can  lead  to  serious  errors. 
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CONCLUDING  REMARKS 


Calculations  have  been  reported  of  the  dissociative  attachment  of  low  energy  electrons  to 
molecular  hydrogen  in  several  vibrational  levels.  The  calculations  are  based  on  fully  ab  initio 
calculations,  obtained  from  the  literature,  of  the  potential  curves  and  electronic  structure  matrix 
elements  that  characterize  the  fixed-R  electron-H2  scattering.  The  present  calculations  have  fully 
included  the  energy-dependent  and  nonlocal  width  and  shift  operators. 

By  monitoring  the  accuracy  of  the  calculations  as  the  shift  and  width  operators  are  reduced 
to  local  form,  we  have  demonstrated  that  a  poor  choice  of  the  local  shift  can  seriously  degrade  the 
accuracy  of  the  results,  even  if  the  correct,  energy-dependent  and  nonlocal  width  is  used. 
Unfortunately,  the  best  choice  of  a  local  resonance  potential  curve,  which  follows  from  the 
optimum  local  shift,  is  not  the  one  available  from  fixed-nuclei  calculations  of  the  resonant  potential 
curve.  Furthermore,  our  results  suggest  that  a  satisfactory  local  approximation  to  the  shift  operator 
may  not  always  exist.  Our  experience  with  H2  suggests  that  an  important  requirement  is  that 
several  vibrational  channels  be  open  at  the  threshold  for  DA.  For  H2, 10  vibrational  channels  are 
available  at  this  energy,  so  the  requirement  is  satisfied.  The  situation  could  be  quite  different  for 
other  systems. 

This  analysis  of  the  results  underscores  the  need  for  detailed  electronic  structure 
calculations  of  the  fixed-R  electron-H2  scattering.  In  particular,  the  energy  dependent  matrix 
element  V(R,e)  coupling  the  resonant  state  with  the  continuum  must  be  obtained  for  a  definitive 
calculation.  The  present  work  has  shown  that  once  the  necessary  fixed-nuclei  information  has 
been  determined,  the  dynamics  equations  involving  energy-dependent  and  nonlocal  operators  can 
be  efficiently  formulated  and  solved. 
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FIGURE  CAPTIONS 


1 .  Comparision  of  the  numerical  values  r(R,e)  =  27tlV(R,e)l2  calculated  by  Miindel  et  al.  with 
the  analytic  fit  used  in  the  present  work. 

2.  Determination  of  the  least  squares  fit  to  A(R).  The  values  of  r(R,e)  =  2jcIV(R,e)I2  shown  in 
Fig.  1  and  calculated  in  Ref.  8  can  be  inverted  using  Eq.(14)  to  give  a  value  of  A(R)  for  each 
(R,e)  pair.  The  range  of  values  of  A  thereby  calculated  at  each  R  is  indicated  by  the  boxes. 

3 .  The  reduced  function  used  to  determine  the  principal  part  integral  defined  by  Eq.(  1 8). 

4.  The  relevant  potential  curves  for  the  dissociative  attachment  of  electrons  to  H2. 

5 .  Calculated  cross  section  for  dissociative  attachment  of  electrons  with  H2(v,J  =  0)  through 
the  2Zu  resonant  state.  The  full  energy  dependent  and  nonlocal  shift  and  width  operators 
were  used  for  these  calculations. 

6.  The  calculated  cross  section  for  dissociative  attachment  of  electrons  with  H2(v  =  0J  =  0)  is 
compared  with  several  previous  calculations  and  with  the  experimental  data  of  Schulz  and 
Asundi. 

7 .  The  peak  cross  section  for  dissociative  attachment  as  a  function  of  the  internal  energy  of  the 
molecule.  The  scale  on  the  top  of  the  figure  shows  the  corresponding  vibrational  quantum 
number  for  J  =  0. 

8 .  Comparision  of  rate  constants  for  dissociative  attachment  with  the  corresponding  calculations 
of  Wadehra  based  on  the  local  complex  potential  model. 
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9.  Comparision  of  the  cross  section  for  dissociative  attachment  calculated  using  the  full  energy- 
dependent,  nonlocal  shift  operator  with  a  calculation  using  the  local  shift  defined  by  Eq.(31). 
The  energy-dependent,  nonlocal  width  operator  was  used  for  both  calculations. 

10.  Comparison  of  two  forms  of  the  local  potential  used  for  the  results  presented  in  Fig.  9.  The 
local  potentials  Wjts  and  Vjocai  are  obtained  by  adding  to  the  potential  Vd  different  local 
approximations  for  the  shift  operator. 
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Figure  1 
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Figure  4 
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